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Abstract 

We review the physics of nuclear matter at high energy density and the experimental 
search for the Quark-Gluon Plasma at the Relativistic Heavy Ion Collider (RHIC). The 
data obtained in the first three years of the RHIC physics program provide several lines 
of evidence that a novel state of matter has been created in the most violent, head-on 
collisions of Au nuclei at y/s = 200 GeV. Jet quenching and global measurements show 
that the initial energy density of the strongly interacting medium generated in the colli- 
sion is about two orders of magnitude larger than that of cold nuclear matter, well above 
the critical density for the deconfinement phase transition predicted by lattice QCD. The 
observed collective flow patterns imply that the system thermalizes early in its evolution, 
with the dynamics of its expansion consistent with ideal hydrodynamic flow based on a 
Quark-Gluon Plasma equation of state. 
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1 Introduction 

Quantum Chromodynamics (QCD) is the fundamental theory underlying the strong interaction 
between quarks and gluons that dictates the structure of hadrons and nuclei. The non-Abelian 
nature of the QCD gauge symmetry gives rise to confinement of quarks and gluons, making 
hadrons the only stable vacuum excitations of the strong interaction. Under extreme conditions 
of high temperature or baryon density the structure of the vacuum is predicted to change, 
undergoing a phase transition that restores the broken symmetries [IJ and frees quarks and 
gluons from confinement [2l!31- Numerical calculation of the equation of state (EOS) in lattice 
QCD theory indeed finds a nearly first-order phase transition at a critical temperature of about 
170 MeV jlj. Such a form of matter, called the Quark-Gluon Plasma (QGP) |Sl!nilZ|, existed in 
the early universe a few microseconds after the Big Bang. Quark matter may still exist today 
in the core of neutron stars, where the baryon density exceeds the critical value of the phase 
transition jH]. 

Soon after the QGP was identified as the excited phase of the QCD vacuum, it was realized 
that high energy collisions of heavy nuclei could create large volumes of matter at high energy 
density [DJ and that the energy density achievable at current or foreseeable nuclear accelerators 
might be sufficient to create the QGP in controlled laboratory experiments. Searches for the 
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QGP using heavy-ion collisions have been carried out over the past three decades at successively 
higher energy facilities, beginning at the Berkeley Bevalac in the early 1980s JU] and continuing 
at the Brookhaven AGS and CERN SPS in the late 1980s ^T]. The Relativistic Heavy Ion 
Collider (RHIC), the subject of this review, was commissioned in 2000. Future experiments are 
planned at even higher energies at the CERN LHC |T2]. 

Experimental results from the AGS and SPS have revealed a rich set of new phenomena 
indicating the formation of dense matter notably the anomalous suppression of char- 

monium jTH] and the broadening or mass-shift of vector mesons |T^. However, it has proven 
difficult to disentangle hadronic contributions to the observed signals, and no clear consensus 
has emerged that a long-lived Quark Gluon Plasma has been observed in the fixed target heavy 
ion experiments at those facilities. 

Experiments at RHIC, the world's first heavy ion collider, have initiated a new era in the 
study of QCD matter under extreme conditions. Nuclear collisions at the highest RHIC energy 
(nucleon-nucleon CM energy y/s^=200 GeV) not only produce matter at the highest energy 
density ever achieved, but also provide a number of rare observables that have not been accessible 
previously. These observables are especially clean probes of the hot and dense matter generated 
in the collision. The wealth of data collected and analyzed in the first three years of RHIC 
operation indicates that a dense, equilibrated system is generated in the most violent, head- 
on collisions of heavy nuclei. The initial energy density probed by hard processes and global 
measurements is estimated to be about two orders of magnitude larger than that of cold nuclear 
matter, well above the critical density for the deconfinement transition predicted by lattice QCD 
calculation. The collective behavior of the observed final-state hadrons provides evidence of early 
thermalization and hydrodynamic flow consistent with a Quark-Gluon Plasma equation of state. 
The matter is evidently a near-ideal, strongly coupled fluid, and is quite different from the ideal 
non-interacting gas expected from QCD at asymptotically high temperature. 

The flrst set of experimental runs to survey the RHIC physics landscape has now been 
completed and it is opportune to pause and examine where things stand. In this review we 
will outline the RHIC heavy ion scientiflc program and discuss the main experimental results 
that support the foregoing conclusions. A large array of data is available, with over 50 Physical 
Review Letters published by the four RHIC experiments thus far. In limited space we cannot 
review aU of the physics topics under study at RHIC. We therefore concentrate on those topics 
for which the data are mature and most clearly address the issues of QGP formation. There 
has also been considerable recent interest in universal properties of high density QCD at low 
Bjorken xbj PZ] and we will discuss RHIC data relevant to this physics. There are a number 
of signiflcant topics that we touch only lightly or omit altogether, among them fluctuations and 
heavy quark and quarkonium production. These areas are developing rapidly and a review at 
this time would soon become outdated. Detailed discussions of all aspects of RHIC heavy ion 
physics can be found in the proceedings of recent Quark Matter conferences [THl ITU] . 

The review is organized as follows: section 2 presents a theoretical overview of the QCD phase 
transitions, together with theoretical considerations for analyzing heavy-ion nuclear collisions; 
section 3 presents the RHIC collider and experiments; section 4 discusses bulk particle production 
and the constraints it places on the collision dynamics; section 5 presents evidence from collective 
phenomena that equilibrium is achieved early in the flreball evolution; and section 7 discusses 
hard probes, concentrating on the theory and measurements of jet quenching. Section 8 presents 
a summary and outlook. 
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2 QCD and the Quark Gluon Plasma 



Within the Standard Model, the strong interaction between quarks is mediated through non- 
Abehan gauge gluon fields, as described by a theory called Quantum Chromodynamics (QCD), 
with 

^QCD = E ^^M^d^ -9A^a^- m,)V^, - - E K^'Fa,,.. (1) 

where Aa's are Gell-Mann matrices in the fundamental representation of S\J (3) and F^^ = d'^A'^— 
A^^ + igfabcA'i^A'^ are gluon field-strength tensors. The seemingly simple theory is very similar 
in form to Quantum Electrodynamics (QED) but it possesses much richer structure because of 
its many symmetries, including the SU{3) gauge symmetry, the approximate chiral symmetry 
for the light quarks that is spontaneously broken in the QCD vacuum. Ua{^) symmetry and 
scale invariance are both broken through quantum interactions, leading to Ua{^) (or chiral) and 
scale (or trace) anomalies. These symmetries and their breaking dictate the structure of the 
vacuum and properties of strongly interacting matter, including the different phases that have 
been the focus of many theoretical and experimental studies since QCD was established as the 
theory of strong interactions. 



2.1 Deconfinement 

Many remarkable features of QCD can be traced to the underlying SU (3) gauge symmetry of the 
strong interaction between quarks and gluons, which both carry color charges. Because of the 
non-Abelian self-interaction among gluons, the color charges at short distance are anti-screened 
due to color diffusion via gluon radiation, leading to a weakening of the coupling constant. This 
asymptotic freedom of the strong interaction |2IIl EI] opens the door for perturbative studies 
of the strong interaction within QCD, including renormalization. The scale dependence of the 
strong coupling, which at the leading order is given by 

4:TT 

"^^^'^ = (ll-|n,)logQ7A|cD' 

has been successfully tested by many experiments |22j. While the small value of the strong 
coupling constant implies weak interaction and thus makes it possible to calculate physical 
processes at short distance, the divergence at small energy scale or long distance signals strong 
interactions and therefore color confinement. Indeed, non-perturbative calculations of the heavy 
quark potential in lattice QCD show a linear potential at long distance with a string tension k ^ 1 
GeV/fm I24j . indicating confinement of quarks to the domain of hadrons in normal vacuum. 
Since the average inter-particle distance becomes smaller at higher temperature or density, the 
interaction among quarks and gluons becomes weaker and confinement will eventually disappear, 
leading to a new phase of matter called the Quark Gluon Plasma (QGP) which is distinctly 
different from the hadronic phase. The energy density of such non-interacting QGP is 



77r2 n^' 

Guf h 16 — 

^ 120 30 



(3) 



with the pressure P = e/3. In contrast, e^r = 3'k'^T'^/30 for a massless pion gas. 

Related to confinement, the renormalization of the coupling constant also breaks the scale 
invariance of QCD, leading to the trace anomaly T'^^ = {as/l27i)F'^'^F^i, + mipip of the energy- 
momentum tensor. The non-vanishing value of the vacuum expectation value of the gluon 
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Figure 1: Left: The energy density as a function of temperature from lattice QCD J^. Arrows 
show the ideal gas values esB from Eq. Q. Right: Deviation from ideal gas EOS (e — 3P)/T'^ 
at fiB=0, 210, 410 and 530 MeV (bottom to top) as a function of T/T^ gZj. 

condensate j2Sl, 

48B = {^F^^F^,) ^ 0.02GeV' (4) 

TT 

implies positive pressure and energy density in the vacuum that confines the weakly interacting 
quarks and gluons inside a hadron according to the MIT bag model [26 . This concept can be 
extended to the equation of state of a non-interacting QGP and massless pion gas, generating 
a first-order phase transition at a temperature Tc ~ 0.725^/^ for a baryon-free system. Such a 
phase transition has indeed been found in numerical simulations of QCD on the lattice. Fig. Q 
shows results from a recent lattice QCD calculation of the energy density as a function of 
temperature jl]. The sharp rise of the energy density at ~ 170 MeV for two light quark 
flavors signals an increase of the effective number of degrees of freedom. The transition becomes 
flrst order for three light flavors, but the order of the phase transition is still not clear for a 
realistic value of the strange quark mass. However, a sharp cross-over is clearly present. 

Shown as horizontal arrows in Fig. ^ (left panel) are the energy densities of a non-interacting 
and massless parton gas as given by Eq. Q. Lattice QCD results clearly deviate from these 
limits, indicating strong interaction among partons even above the phase transition temperature 
Tc. Higher order corrections to the perturbative calculation in the weak coupling region do not 
improve agreement with the lattice results. For reasonably small values of coupling constant g or 
Os, the perturbative expansion series up to order g'^ is found not to converge |2S]. This indicates 
that massless partons are not a good basis for perturbative expansion because of the contribution 
of soft modes with momentum k ~ gT. Resummation of hard thermal loops (HTL)^^ leads 
to a picture of quasi-particle modes dominating the interaction in an interacting QGP. Working 
on the basis of quasi-particles and solving Dyson's equation in a self-consistent resummation, 
an equation of state (EOS) is obtained that is very close to that of lattice QCD results at 
temperatures above but close to the phase transition temperature [5U] . 

Since the EOS of an ideal massless QGP is P = e/3, a measure of deviation from an ideal 
gas EOS is e — 3P. Lattice QCD results [^Z] (Fig. [TJ right panel) indicate the presence of 
strong interactions around the phase transition temperature. Since the sound velocity is given 
by = dP/de, the softening of the EOS around the phase transition region would slow down 
hydrodynamic expansion, with signiflcant consequences for the evolution of the dense matter 
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and the relation of observed transverse energy production to the initial energy density prior to 
the start of hydrodynamic expansion. 

In addition to large energy density, other properties of the quark gluon plasma will be quite 
different from those of a hadron gas or normal nuclear matter. Hadrons as bound states of 
quarks and anti-quarks do not exist in the plasma. Lattice studies of the heavy quark potential 
find that the confining linear potential in vacuum disappears once the temperature is above the 
phase transition temperature and what is left over is a screened Coulomb potential jHl]. This 
is why heavy quarkonium was proposed as a good probe of the quark gluon plasma formed in 
high-energy heavy-ion collisions since the quarkonium bound state will be dissolved into 
an unbound quark- ant iquark pair once the binding potential is screened in the plasma. Though 
various theoretical estimates also predict quarkonium break-up through hadronic interactions 
|331 ES] , thermal gluon exchange screening of potential in a quark gluon plasma remains the 
most efficient way to suppress quarkonium. 

Simulations on the lattice are difficult for a system at finite baryon chemical potential ij,b- 
However, significant progress in this area has been made in recent years. The phase transition 
is found to change from first-order at high baryon density to a cross over at lower values of 

EZj , though the position of the end-point is still uncertain and depends on the values of 
quark masses. Such an end-point in the phase diagram of QCD matter could have important 
effects, for instance increased baryon number fiuctuations in heavy-ion collisions at low energy 
if the end-point is crossed over in the evolution of the fireball. 

At high baryon densities and zero or low temperatures, the attractive interaction between 
quarks in the color-anti-triplet and spin- zero channels could lead to pairing of quarks in iso-singlet 
spin-zero and color-anti-triplet states on the Fermi surface, thus forming diquark condensates in 
two fiavor quark matter jnilESI- With three fiavors of quarks, the Cooper quark pairs cannot 
be iso-singlet. The attractive interaction is found to be favored in the channel where colors and 
fiavors are locked [10]. Such color-fiavor locking leads to a plethora of new structures in the 
phases of dense quark matter which may generate novel effects in neutron stars where such high 
densities could be reached. 

2.2 Chiral Symmetry Restoration 

In addition to the SU{3) local gauge symmetry, the QCD interaction as given by Eq. has 
the global symmetries of SU{3) x ^[/^(S) x U{1) x Ua{1) in the quark sector. The global 
U{1) symmetry is responsible for baryon number conservation while chiral symmetry, which is 
approximate for vanishing quark masses, also leads to many unique properties of QCD. 

If the quark fields are decomposed into left and right chirality components, ipL,R = (1 ^75)"^) 
the QCD Lagrangian with three massless quarks is invariant under the transformation. 



According to the Noether theorem, this will generate two kinds of conserved currents or their 
combination, vector and axial-vector currents. 



Identifying each current with a physical state having the corresponding quantum number results 
in degeneracy between scalar and pseudoscalar as well as between vector and psudovector mesons. 



(5) 




(6) 
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The absence of parity doublets led to the the discovery of the spontaneous breaking of chiral 
symmetry by the QCD vacuum, manifested by the non-vanishing of the quark condensate 

= {i,Ri,L + ^l^k) = -(240MeV)=^ , (7) 

which is directly related to /^r, the pion decay constant. One of the consequences of this spon- 
taneous symmetry breaking is the existence of massless Goldstone bosons. The absence of the 
ninth Goldstone boson is due to the breaking of Ua{^) symmetry by quantum correction or an 
anomaly, 

d,A>', = -ia.FrF,,,., (8) 

iOTT 

where Fa^^u = e^yaisF^^ . The non-vanishing vacuum expectation value of this anomaly is related 
to what is known as the topological susceptibility of the vacuum. 

All of these order parameters, the quark and gluon condensates and the topological suscep- 
tibility, are believed to be related to the gluonic structure of the vacuum |1T]. They dictate all 
hadronic properties - their mass spectra, decay widths and constants. The medium modifica- 
tion of these condensates at low temperature and density is calculable within the framework of 
chiral perturbation theory 02] • At higher temperature and in the quark gluon plasma phase, 
they are expected to vanish. The quark condensate in lattice QCD calculations |1| is shown 
to disappear above the QCD phase transition temperature characterized by a sharp cross-over, 
which coincides with the cross-over of the expectation value of the Polyakov loop, a measure of 
deconfinement. At these temperatures and densities, chiral symmetry is completely restored and 
[/^(l) symmetry is also partially restored. Hadronic properties will be very different from those 
in vacuum. There will be mixing of vector and axial vector currents and all the chiral multiplets 
will become degenerate. 

There are many consequences of the restoration of chiral and Ua{^) symmetry in the dense 
medium. The vanishing value of the topological susceptibility gives rise to a reduction of the 
ninth Goldstone boson, leading to possible enhancement of rj mesons |431 144j . A rapid cooling 
of the dense medium initially in a chirally symmetric state could introduce instabilities into 
the evolution of chiral fields leading to the amplification of soft modes j33], though a realistic 
treatment of the expansion and cooling of a finite system can prevent the creation of such 
disoriented chiral condensates (DCC) in large domains [IHl- Perhaps the most promising signals 
of chiral symmetry restoration are modifications of hadron properties in medium, i. e. the masses 
and widths of vector mesons and axial vector mesons via mixing, as manifested in the soft 
dilepton spectra |17j. The enhanced soft dilepton yield below the p mass region seen by the 
CERES experiment at SPS may be attributable to the medium modification of the p meson. 
During the QGP phase, annihilation of thermal quarks and antiquarks could also contribute to 
the underlying dilepton spectra in the intermediate mass region. 

2.3 Perturbative QCD 

The most important practical consequence of asymptotic freedom for the SU (3) gauge interaction 
in QCD is the success of perturbative QCD (pQCD) in describing various processes involving 
strong interactions at short distance. Because of the small coupling constant, it is possible to 
make a systematic expansion of physical cross sections in the strong coupling constant when 
the momentum transfer involved is large. Because of color confinement in vacuum, both the 
initial and final observed particles will be hadrons and therefore will involve strong interaction 
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at long distance, which is not calculable within the framework of the perturbative expansion. 
However, it has been proven to at least leading power correction (l/Q^) that the cross section can 
be factorized into short-distance parts calculable in pQCD and non-perturbative long distance 

parts iiHiiniEni- 

The long distance parts normally involve hadronic wavef unctions, parton distributions, and 
hadronization. Though they cannot be calculated perturbatively, their matrix element definitions 
are universal, independent of specific processes. If measured in one process they can be applied 
to another process; therein lies the predictive power of pQCD. The renormalization of these 
non-perturbative matrix elements due to initial and final state radiation in pQCD leads to the 
Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) evolution equations [B^ESESI- Given the 
experimental measurements of these distributions at one scale, the DGLAP equations predict 
their evolution to higher momentum scales. For example, final-state radiation leads to the 
DGLAP evolution equations for parton fragmentation functions Da->h{zh, Q"^), 



ainQ2 27r z 



'-fq^qg{z)Dq^h{Zh/z, Q^) + -fq^gq{z) Dg^h{Zh/ Z , Q^) 
■2nf 

J2 lg^qqiz)Dq^h{Zh/ Z, Q"^) + 1 g-,gg{z) D g^h{Zh / Z, Q^) 
q=l 



(9) 



where •ya-^bciu) are the splitting functions of the corresponding radiative processes j^ni- These 
evolution equations have been tested extensively against experiment and can now even be used 
to measure the scale-dependence of the running strong coupling constant |22]- 

The perturbative QCD parton model is based on this factorization picture of hard processes. 
In this model the cross section of a typical hard process can be expressed as the convolution of 
initial parton distributions, perturbative parton scattering cross sections and the parton frag- 
mentation function. Since hard processes happen on a short time scale in the very earliest stage 
of high-energy heavy-ion collisions, they probe of the bulk matter that is formed shortly after 
the collision. The pQCD parton model serves as a reliable and tested framework for the study of 
these hard probes. Proposed hard probes in high-energy heavy-ion collisions include Drell-Yan 
dileptons from quark- antiquark annihilation, direct photons from quark and gluon Compton 
scattering, heavy quark production, and high pj- jets from hard parton-parton scattering. In 
this review, we will focus on the physics of jet propagation in the dense medium. 



3 The Relativistic Heavy Ion Collider 

The Relativistic Heavy Ion Collider (RHIC) jSl] at Brookhaven National Laboratory is the 
world's highest energy accelerator of heavy nuclei and the world's first polarized proton collider. 
In this section we describe the accelerator and the experiments, together with some theoretical 
considerations important for the analysis of ultrarelativistic heavy ion collisions. 

3.1 The Accelerator 

The collider consists of two independent, concentric acceleration and storage rings, with a cir- 
cumference of 3.8 km. All storage ring magnets are superconducting, cooled to 4.2 K by a 25 kW 
helium refrigerator. There are six intersection points, of which four are currently instrumented 
with experiments. RHIC can store and collide beams with masses ranging from protons to the 
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Au+Au 


P+P 


Beam energy 
Mean Luminosity at top energy 
Bunches per ring 
Luminosity lifetime 

13* 


30 — > 100 GeV per nucleon 
2 X 10^6 cm-^s-i 

60 120 
~ 10 hours 
10 ^ 1 m 


30 ^ 250 GeV 
1.4 X 10^1 cm-^s-i 
60 ^ 120 
>10 hours 
10 ^ 1 m 



Table 1: Main design specifications of RHIC for Au+Au and p+p collisions |54j . 



heaviest stable nuclei. Due to the independence of the rings RHIC has great flexibility to collide 
beams of unequal masses, such as protons or light ions with Au ions. The top collision energy 
for the heaviest nuclear beams is y^s^=200 GeV per nucleon pair, while the top energy for p + p 
is yi=500 GeV. 

The layout of RHIC is shown in Fig. |21 Heavy ion beams originate in a pulsed sputter 
source and are accelerated successively by a Tandem van der Graaf accelerator. Booster Syn- 
chrotron, and the Alternating Gradient Synchrotron (AGS), where they are accelerated to 10.8 
GeV/nucleon, fully stripped of their electrons, and injected into RHIC. Polarized protons orig- 
inate in a 200 MeV Linac and are accelerated by the Booster and the AGS to 24.3 GeV for 
injection into RHIC. Polarization of protons is maintained by use of Siberian Snakes [55j. The 
physics of the RHIC Spin program is beyond the scope of this review; a recent review and status 
report can be found at [SHIEZI- 

Acceleration and storage in RHIC utilizes two Radio Frequency (RF) systems, one at 28 MHz 
to capture the AGS bunches and accelerate to top energy, the other at 197 MHz to provide a short 
collision diamond (cr^ ~ 25 cm) for efficient utilization of the luminosity by the experiments. 
The synchotron phase transition of the RHIC lattice is at 7t = 24.7, meaning that all ions 
except protons pass through the beam instability at transition. 

The main performance specifications of RHIC are given in Tabled For light ions (yl<100), 
the luminosity is limited by beam-beam hadronic interactions, whereas for heavier ions the 
luminosity lifetime is limited by intrabeam (intra-bunch) scattering [58] • Other processes which 
significantly limit the luminosity for heavier ions result from their intense electromagnetic fields 
at high energy: Coulomb dissociation and spontaneous electron-positron production followed by 
electron capture |2H] • The luminosity for beams of heavy nuclei scales with beam energy as 7^ . 

Table 121 lists the beams, energies and integrated luminosity for the RHIC runs to date. The 
recently completed runs labelled "2004" have not generated scientific publications as of this 
writing. While the luminosities accumulated for heavy ion running appear at first sight to be 
small, nuclear geometric effects amplify hard cross sections approximately as the product of the 
nuclear masses AB (Sect. 13. 3^ . The rightmost column converts the integrated luminosity for 
heavy ions to the equivalent for p + p collisions that would generate the same yield for a process 
whose cross section scales as AB. In the later runs the store-averaged luminosities achieved for 
Au + Au collisions exceed the RHIC design specifications in Tabled 

3.2 The Experiments 

The RHIC beams are brought into head-on collision at the intersection regions. The final dipoles 
of the lattice are approximately ±10 m from the collision diamond, with the intervening space 
free for detectors. Currently four of the six intersection regions are instrumented, with two major 
detectors (PHENIX, STAR) and two smaller ones (BRAHMS, PHOBOS). 
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Figure 2: The RHIC accelerator complex. 



Beams 


^/^(GeV) 


/m(nb-i) 


p + p equivalent / Cdt (pb ^) 


Au+Au 


20 


small 


small 


Au+Au 


130 


0.02 


0.8 


Au+Au 


200 


0.24 


7.9 


p + p 


200 


1600 


1.6 


d+Au 


200 


75 


3.0 


Au+Au (2004) 


200 


~2 


80 


Au+Au (2004) 


62.4 


-0.05 


2 



Table 2: Beam and energy combinations run at RHIC to date. / Cdt denotes integrated luminos- 
ity delivered to all four experiments, "p + p equivalent" is the luminosity scaled by the number 
of binary collisions (note change of units between columns 3 and 4). Runs labelled "2004" were 
recently completed and have not generated physics publications as of this writing. 
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All four experiments contain identical Zero Degree Calorimeters (ZDC) [59j, which are used 
for triggering, luminosity monitoring, and event characterization. The ZDCs are compact tung- 
sten/fiber hadronic calorimeters situated immediately downstream of the machine dipole mag- 
nets that define the interaction region, with acceptance 2.5 mrad centered on the beam direction. 
The energy flux into this acceptance is dominated in heavy ion collisions by non-interacting spec- 
tator neutrons. The ZDC energy resolution is sufficient to discriminate individual beam-velocity 
neutrons [22] • The luminosity for Au + Au collisions has been measured with 10% precision 
using a Vernier scan of ZDC coincidence rates [HIH Ifilj. 

PHENIX The PHENIX experiment is designed to make precision measurements of a wide 
variety of observables, sensitive to multiple time scales in the evolution of heavy ion collisions. 
Special emphasis is put on rare signals (direct photons, lepton pairs, J/ip and T families, jet 
fragments) that probe the system at the earliest, hot and dense phase. Measurement of lepton 
pairs at low px is a promising tool for studies of chiral symmetry restoration. Inclusive measure- 
ments and correlations of identified hadrons at high pt are sensitive to partonic interactions in 
the medium, and at lower px probe the late, hadronic gas stage of the collision. 

The PHENIX detector jH21; shown in Fig. El consists of four independent spectrometers, 
two at midrapidity for charged hadrons, electrons, and photons, and two at forward rapidities 
for muons. Each has an acceptance of ~ 1 steradian. The midrapidity spectrometers have 
an axial magnetic field, with tracking for momentum measurements supplied by drift and pad 
chambers. Charged particle identification over a broad momentum range is provided by Time of 
Flight (TOF), Ring Imaging Cerenkov (RICH), and Time Expansion Chamber (TEC) detectors, 
giving proton identification up to pt=5 GeV/c. Electrons and photons are measured in highly 
granular lead-scintillator and lead glass calorimeters (EMC). The combination of EMC, RICH 
and TEC provides a hadron background rejection factor for electron measurements of 10^ over a 
wide momentum range. The forward muon spectrometers have acceptance for J/ip ^ fi~^fi~ of 
—2.25<y< — 1.15 and 1.15<|/<2.44. The muon arms have a radial magnetic field, with tracking 
based on drift chambers followed by a muon identifier consisting of alternating layers of steel 
absorber and streamer tubes. Pion contamination of identified muons is ~ 3 x 10^^. 

STAR The STAR experiment has broad physics reach, covering a wide variety of hadronic and 
leptonic observables. Large acceptance enables measurement of a large fraction of the thousands 
of charged hadrons produced in a heavy ion collision, to measure correlations and to search for 
rare or subtle non-statistical fluctuations indicating new physics. STAR has robust capabilities 
to trigger and measure high p^ observables (hadron yields and correlations, electrons, photons, 
and jets) to investigate partonic interactions in dense matter. Measurements of J/ip and T will 
probe deconfinement directly. 

The STAR detector [HSI is shown in fig. 0] STAR is based on a warm solenoidal magnet, with 
radius 260 cm and maximum field strength 0.5 Tesla, surrounding a variety of detector systems. 
The main tracking device is a solenoidal Time Projection Chamber (TPC) with radius 200 cm 
and full azimuthal acceptance over |?7|<1.4. Additional tracking is provided by inner silicon drift 
detectors at midrapidity and forward TPCs at 2.5<|r7|<4. Photons and electrons are measured 
by Barrel and Endcap Electromagnetic Calorimeters (EMC), with full azimuthal acceptance 
over — 1.0<r7<2.0. Particle identification is carried out using specific ionization (dE/dx) in the 
TPC gas, time of flight, reconstruction of displaced vertices for weakly decaying particles, and 
combinatorial invariant mass methods. The identiflcation of strange baryons and mesons has 
been made up to Pt=Q GeV/c and of charmed mesons to 10 GeV/c, with the measurements 
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currently limited by statistics. Fast triggering utilizes the ZDCs, forward scintillators (Beam- 
Beam counters), a barrel of scintillator slats surrounding the TPC, and the EMC. 

BRAHMS The emphasis of the BRAHMS experiment is on high precision inclusive measure- 
ments and small-angle correlations of primary hadrons over very broad phase space. BRAHMS 
indeed has very large coverage, extending for charged pions to rapidity ?/ = 4 [Au-heam rapidity 
is 5.37). The BRAHMS detector jMj, shown in fig. consists of two independent charged 
particle spectrometers: the Forward Spectrometer (FS), with acceptance 0.8 msr measuring mo- 
menta pr<35 GeV/c for angles relative to the beam 2.3<6'<30 degrees, and the Mid- Rapidity 
Spectrometer (MSR), with acceptance 6.5 msr and angular coverage 30<6'<95 degrees. Mo- 
mentum measurements and particle identification utilize Time Projection Chambers (TPCs) in 
conjunction with Time of Flight hodoscopes and Threshold and Ring Imaging Cerenkov detec- 
tors. Triggering and global even characterization of heavy ion events are carried out using a 
mid-rapidity multiplicity array, forward scintillator detectors, and the ZDCs. 

PHOBOS PHOBOS is designed to carry out a very general search for a priori unknown and 
potentially rare signatures of new physics, requiring a very large acceptance device that detects 
almost all charged particles in each event for a large fraction of all inelastic collisions. Virtually 
the entire RHIC phase space is covered by the PHOBOS multiplicity measurement, with the 
high trigger and recording rate allowing offline searches for unusual events or rare fluctuations. 
Large-scale phenomena in heavy ion collisions may generate effects at very low pt, requiring the 
spectrometers to measure down to pt ~30 MeV/c. The PHOBOS spectrometers also measure 
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at sufficiently high px to be sensitive to jet-related observables. 

The PHOBOS experiment shown in fig. IHl is based almost entirely on silicon pad detec- 
tors. It consists of a multiplicity array covering 11 units of pseudorapidity, a finely segmented 
vertex detector, two small-acceptance midrapidity spectrometers, and trigger detectors. The 
spectrometer arms utilize a warm dipole magnet of strength 1.5 Tesla-m. Particle identification 
is based on measurements of time of flight and energy loss in the silicon. Special emphasis is put 
on measurements at very low transverse momentum, requiring a thin beam pipe and minimal 
material in front of the first tracking planes. 




Figure 6: The PHOBOS detector. 



3.3 Geometric Aspects of High Energy Nuclear Collisions 

In proton-nucleus and nucleus-nucleus collisions, multiple scattering occurs at both the hadronic 
and the partonic level. Multiple scattering influences many aspects of the dynamics of high 
energy nuclear collisions, in particular the initial parton scattering responsible for bulk particle 
production and the spectrum of rare hard processes. It depends on the geometry of the nucleus- 
nucleus collision, which also dictates the geometry of the produced dense matter. Determination 
of the collision geometry is a key element in the study of heavy-ion collisions. 

Multiple Scattering and the Glauber Model 

Current theoretical treatment of multiple scattering is based on the Glauber model jnHl- In this 
model, a hadron-nucleus collision is considered to be a series of multiple hadron-nucleon scat- 
terings. Neglecting the difference between the hadron and its excited states between successive 
scatterings and utilizing the forward peak of high energy hadron-nucleon elastic scattering, the 
total cross section of h + A collisions is [nH] 




where tA{b) = TA{b)/A and is the nuclear thickness function defined as a line integral over 

the nuclear density pA, 

TA{h) = / dzpA{b,z). (11) 

J —oo 
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The exponential factor in the above cross section comes from the interference between different 
scatterings, assuming a Gaussian form of the nuclear density distribution. The longitudinal 
momentum transfer is related to the transverse momentum transfer qrpi of each scattering, 



(12) 



and determines the coherence length 1^ = ^/^Qzi- For high energy scattering with small trans- 
verse momentum transfer, the coherence length is much longer than the nuclear size Ra and 
hadron-nucleus collisions become coherent. The hadron-nucleus cross section is then given by 
the familiar Glauber formula. 



cThA = I d''b{l-[l-ahNtAib)]''} . 



(13) 



For large nuclei A ^ 1, the integrand can be approximated by an exponential 1— exp[— (T/iArTA(6)]. 
For hadronic scattering with large cross section ahN, this cross section is mainly determined by 
the nuclear density distribution and is denoted the geometric cross section. This approach is 
most relevant for soft processes. 

For hard processes with large transverse momentum transfer and small cross section, the 
coherence length becomes much smaller than the intra-nucleon distance in the nucleus. In this 
case all interference terms drop out and the cross section results from the incoherent superposi- 
tion of nucleon-nucleon collisions. The h + A cross section is then directly proportional to that 
foTh + N collisions: 



hard 
hA 



hard 
hN 



/I hard 



(14) 



This expression is proportional to the thickness function, which represents the number of hadron- 
nucleon collisions in a h + A collision. As shown in Fig. [7| (left panel), experimental cross sections 
of Drell-Yan dilepton production with large invariant mass at the CERN SPS indeed scale linearly 
with the atomic number to good precision in p + A collisions. 
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Figure 7: Scaling of hard and soft processes in nuclear collisions at the CERN SPS. Left: Drell- 
Yan production cross section in p+A collisions normalized by atomic mass A (Eq. IT^ . from 
NA50 [67J. Right: charged hadron pseudorapidity distribution in central collisions of equal mass 
nuclei scaled by the number of nucleon participants (~ 2 A), from NA49 [HH]- 
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This formula can be extended to the case oi A + B colhsions by replacing the thickness 
function in h + A collisions with the nuclear overlap function, 

TAB{b) = I dhTA{s)TB{b- s), (15) 
which represents the number of binary nucleon-nucleon collisions per unit cross section, 

dNun/d% = TAB{h), (16) 

m A + B collisions at fixed impact parameter h. 

Since the hard processes are incoherent, their cross sections m. A + B collisions should be pro- 
portional to the number of binary collisions Ni,in. We will see below that certain hard processes 
with large momentum transfer provide sensitive probes of the matter generated in nuclear 
collisions through final-state interaction of their reaction products with the medium. Deviation 
from binary scaling of the cross section therefore indicates novel nuclear effects and provides an 
experimental observable to quantify such effects. Isolation of final-state from initial-state effects 
requires comparison of systems in which final state effects are expected to be present or absent. 
This strategy will play an important role in our discussion of hard processes. 

In contrast to hard processes, soft processes typically have large cross section and coherence 
lengths much larger than the nuclear size. The inclusive cross section and total hadron multiplic- 
ity are then expected to scale as the number of participating ( "wounded" ) nucleons (Wounded 
Nucleon Model [HH])- In the Glauber multiple scattering model, the average number of wounded 
nucleons at a fixed impact-parameter in A + B collisions is 



Npart{h) = J d'sTAis) [l - e""^^^^^^-"^] + J d'sTsis- h) [l - e"'^'^^^-^^^-''^] (17) 

Fig. [3 right panel, shows the charged particle pseudorapidity distribution in central collisions of 
equal mass nuclei measured at the CERN SPS, which indeed is seen to scale for massive nuclei 
as the number of participating nucleons. A particle production model embodying such coherent 
processes is the string model implemented in HIJING Monte Carlo model ^fOll^. In this model, 
a wounded nucleon becomes an excited string, with the number of produced hadrons insensitive 
to the number of scatterings suffered throughout the multiple scattering process. This leads to 
hadron multiplicity from soft interactions proportional to the number of participant nucleons. 

We turn now to application of the Glauber model to data analysis. The impact parameter 
h is not observable and measurements necessarily integrate over a finite interval of h. Event 
geometry is tagged by observables correlated with impact parameter, such as charged particle 
multiplicity, whose distributions are binned into percentiles of the total interaction cross sec- 
tion (Fig. in}- Using the Glauber model, differential cross section distributions da/dNun and 
da/dNpart are calculated and the weighted mean values {Nhin) and (Npart) are found within the 
same percentile bins of total cross section. The bins of the measured distribution and model 
calculation are equated, under the assumptions that the event-tagging observable varies on av- 
erage monotonically (but not necessarily linearly) with impact parameter and that fluctuations 
generate negligible mixing of the bin populations. 

The nuclear density in the calculation is taken to be spherically symmetric, with Woods- 
Saxon radial dependence: 

p^(r) = (18) 
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Typical parameters for Au nuclei are density po=0.169/fm^ and surface thickness a = 0.535 ± 
0.027 fm, the latter derived from electron scattering data[72]. The charge radius Tq = 6.38 fm is 
usually increased to ~6.5 fm to account for the neutron skin thickness. 

Two methods are used in the literature to calcuate Nbin and Npart- the Optical and Monte 
Carlo Glauber approaches (see Appendix A of [721 references therein.) The Optical Glauber 
approach is based on a smooth nuclear matter distribution and numerical evaluation of the 
analytic Glauber integrals. The Monte Carlo approach is based on the random distribution of 
nucleons according to the Woods-Saxon density, with nuclear collisions at a given impact param- 
eter modelled by the incoherent interaction of all nucleon pairs. For central Au+Au collisions at 
RHIC energies the practical difference between these approaches is negligible compared to other 
uncertainties of the measurements, but for peripheral collisions the differences can be significant. 
Assessment of this uncertainty must be made when interpreting data that incorporate Glauber 
calculations. 




Figure 8: Charged particle multiplicity distributions from 200 GeV d + Au collisions, from STAR 
|74j . Measurement is for —3.8<ri< — 2.8 (Au-beam direction). Centrality selection and Glauber 
calculations (histograms) described in text. 



As a gauge of the accuracy of the Glauber model. Figure |S1 compares the results of a Monte 
Carlo Glauber calculation to charged particle multiplicity distributions in 200 GeV d + Au col- 
lisions |73]. The measurements are at forward rapidity (— 3.8<?7< — 2.8, Au-heam direction), 
both for minimum bias d + Au collisions and for peripheral d + Au collisions which have a 
beam-rapidity neutron detected in the deuteron beam direction ("ZDC-d"). The calculated 
multiplicity distributions result from the convolution of the forward charged multiplicity distri- 
bution measured in 200 GeV p + p interactions 75\ with the Npart distribution from the Glauber 
model. Both the minimum bias and peripheral multiplicity distribution measurements are well 
reproduced by the model. In addition, the cross for such peripheral collisions is calculated to 
be (18 ± 3)% of the minimum bias cross section, in good agreement with the measured fraction 
(19.2 ± 1.3)%[71j. Figure IHl demonstrates that the Glauber approach provides a sound basis for 
modelling geometric effects in nuclear collisions at RHIC energies. 
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Centrality tagging in Au + Au collisions 

Fig. ini left panel, shows the measured correlation in 200 GeV Au + Au collisions between ZDC 
energy (spectator neutrons) and charged particle multiplicity within 3.1<|?7|<3.9, from PHENIX. 
In the most peripheral (large impact parameter) collisions, the number of forward neutrons and 
the total multiplicity are both small. Both quantities increase for decreasing impact parameter, 
while for the most central collisions the number of spectator neutrons is again small while the 
multiplicity is large. The correlation between these two geometry-sensitive observables is seen 
to be strong. The figure also illustrates the sorting of events into centrality bins corresponding 
to percentile intervals of the cross section, with 0-5% indicating the most central collisions. 




Figure 9: Event centrality characterization. Left: scaled ZDC energy (vertical) vs. forward 
charged multiplicity (horizontal), from PHENIX [ZE]- Right: distribution of charged multiplicity 
at mid-rapidity, from STAR [77j. Both plots are for minimum bias Au+Au events at y^s^=200 
GeV and show the division into centrality bins corresponding to percentile intervals of the total 
cross section. 

The ZDC energy is small for both the most central and most peripheral collisions, and this 
ambiguity limits its utility as a centrality tag. For many applications the multiplicity distribution 
alone suffices as a centrality tag, as shown in Fig. right panel. The shape of the distribution 
is dominated by the nuclear geometry, with the tail at the highest multiplicity governed by 
multiplicity fluctuations within the finite measurement aperture for the most central collisions. 

We conclude this section with a discussion of the impact parameter dependence of Nun and 
Npart- Figure HUl left panel, shows this dependence for Au + Au collisions at RHIC energies, 
indicating the strong bias towards central collisions of binary-scaling processes. For an A^^bm- 
scaling process, we define the fraction of the its total cross section contained in events with 
impact parameter h<hc as [75]: 

fAB = ^ f^hdhTABih). (19) 
AB Jo 

Figure HUl right panel, shows Jab as a function of fgeo, the fraction of the total hadronic in- 
teraction cross section contained in the same impact parameter interval. The A"bj„-scaling cross 
section weights strongly towards central (small impact parameter) collisions, due purely to nu- 
clear geometry. As a rough rule of thumb for binary-scaling processes in symmetric heavy ion 
collisions, 40% of the cross section is contained in the 10% most central events. 
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Impact parameter (b) f 



Figure 10: Left: Nbin and Npart vs. impact parameter for Au+Au at RHIC energies [7^]. Right: 
fraction of binary-scaling cross section vs fraction of hadronic interaction cross section; figure 
from [TH]. Central events correspond to small fgeo- The curves are for (upper to lower) Au+Au, 
Ag+Ag, Cu+Cu, Al+Al, and 0+0 collisions. 

4 Bulk Particle Production and Initial Conditions 

The formation of a Quark-Gluon Plasma in high-energy heavy-ion collisions depends critically 
on the initial parton production at the earliest stage of the reaction. A number of proposed 
signals of the QGP are sensitive to the early conditions and provide direct measurements of the 
initial density and other properties of the dense matter. Alternatively, global observables based 
on bulk particle production, such as the rapidity density of hadron multiplicity and transverse 
energy, provide a coarse-grained view of the early dynamics. These observables may be related to 
the entropy and energy densities early in the collision evolution, and as such provide meaningful 
constraints on the initial conditions that complement the measurements via direct probes. They 
can also provide tests of particle production models and provide a more complete picture of the 
dynamics. Our first discussion will focus on properties of bulk particle production: multiplicity 
and transverse energy pseudorapidity (rj) densities and net baryon number distributions. 

4.1 Particle Production Models 

Prior to RHIC startup, theoretical model estimates of the initial conditions and bulk particle 
production varied over a wide range, due to uncertainties in modeling soft hadron production 
and the interplay between soft and hard processes [201 • These models included a pure pQCD 
parton model a pQCD parton model in combination with string model (TH IH^l, and a 
classical Yang-Mills field model |83|. The uncertainties were significantly reduced with the first 
publication of RHIC data on charged hadron multiphcity [HI]. New calculations based on the 
initial state gluon saturation model [HS] also describe the energy and centrality dependence well. 
Here we provide brief descriptions of three typical models before we present the experimental 
data and discuss their implications. 

Two-component model Mini-jet production in a two-component model was proposed long 
ago to explain the energy dependence of the total cross section |Sni and particle production 
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|87j in high-energy hadron colhsions. This approach was incorporated into the HIJING model 
|7m ITT] to describe initial parton production in high-energy heavy-ion collisions. In this two- 
component model, the nucleon-nucleon cross section at high energy is divided into collisions with 
or without hard or semi-hard jet production processes. The jet cross section ajet is assumed to 
be given by the pQCD parton model. However, the differential jet cross section has an infrared 
divergence as the transverse momentum of the jet goes to zero. An infrared cut-off scale po is 
introduced to separate hard processes from soft processes that are not calculable in pQCD. The 
soft interaction cross section cTgoft is a model parameter. The two parameters, o"soft and po, are 
determined phenomenologically by fitting the experimental data of total p + p{p) cross sections 
within the two-component model [70} [7T|. 

The cut-off scale po separating non-perturbative and pQCD processes could in principle de- 
pend on both energy and nuclear size in A-l-A collisions. Using the Duke-Owens parameterization 
of parton distributions in the nucleon [SH], an energy- independent cut-off scale po = 2 GeV/c, 
and soft cross section cXsoft, the HIJING model can describe well the experimental data on total 
cross sections, hadron multiplicity, and other observables from p + p{p) collisions . The de- 
fault HIJING prediction [90| agreed with the first data on dNch/ dri{\ri\ < 1) of central Au + Au 
collisions at = 130 GeV jHl. 

In applying the two-component model to nuclear collisions, multiple mini-jet production is 
assumed to be incoherent and thus is proportional to the number of binary collisions Nun- The 
soft interaction is however coherent and proportional to the number of participant nucleons Npart, 
according to the Wounded Nucleon Model Assuming no final state effects on multiplicity 
from jet hadronization, the rapidity density of hadron multiplicity in heavy-ion collisions is then 



where cr^fi^s) is the averaged inclusive jet cross section per nucleon-nucleon interaction in AA 
collisions. The average number of participant nucleons and number of binary collisions for a 
given impact parameter can be estimated via Glauber model simulation. Since the parameters 
{n)s and {n)h are determined from p + p{p) collisions, the only uncertainties are due to nuclear 
effects on a^^{s), such as parton shadowing. 

As an alternative to Eq. (j20|) . a simpler parameterization of the "two- component" model is 
given by [HI] 



where is the multiplicity in p + p collisions and x is the fraction of sources scaling as hard 
collisions, related to a^^{s)/ain- 

A final-state parton cascade and hadronic rescattering processes can also be introduced into 
such a pQCD-based model. The AMPT model[H2] uses HIJING for the initial conditions and 
includes parton and hadronic rescattering after-burners. The final-state rescattering is found 
not to change the bulk particle production significantly. 

Final State Saturation A mini-jet pair with transverse momentum po has intrinsic transverse 
area of tt/pq. In A + A collisions at high energy, independent production could result in multiple 
mini-jet production within this area, which is quantum mechanically disallowed. This sets 
the limit on the number of independent mini-jet pairs within a total transverse area 7^R\ as 
N{pq) = pIR\. The pQCD parton model gives the number of scattered partons above transverse 




m 



(20) 




(21) 
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momentum cutoff po as N{pq) = TAA{b)(Tjet{po)- In the final-state saturation model (EKRT JHT), 
the saturation scale Psat is then defined by the self-consistent solution of 



TAA{b)a^et{Psat) = Psat Ra- (22) 

Below the saturation scale, jet production is correlated and the divergent mini-jet cross section 
will be regulated. The EKRT final-state saturation model neglects minijets with pt < Psat and 
assumes that the produced parton density (mostly gluons) and the transverse energy density are 
dominated by minijets above psat- 

Numerical calculation -Sl^ of Eq. ^ yields Psat = 0.21v4°-i3(^yi)0.i9 GeV/c, with initially 
produced gluon multiplicity = 1.38A^'^'^{y/s)^'^^. The calculated ratio of energy density to 
multiplicity density is found to be very similar to that of an ideal gas of bosons, e/n ~ 2. TOT, 
so that for a wide range of A and y/s the gluon gas is generated in a thermalized distribution in 
this model. 

Assuming boost-invariant adiabatic expansion and proportionality between the produced 
parton and the observed hadron multiplicities, the charged hadron density in central collisions 
is given by|HI]: 

^{b = 0) ^ -lA6A''-'^{^f\ (23) 
di] 3 

While this expression only applies for symmetric geometry, it has been extended to non-central 
collisions [OB] and the centrality dependence may be approximated by replacing A with {Npart) 
in the above equation. 

Initial State Saturation More commonly, parton saturation refers to high density effects in 
the initial state 94, .£5^ . For an elementary probe (a virtual photon in deep inelastic scattering, 
a projectile parton in hadronic collisions) interacting with a nucleus of mass A, the coherence 
length of the interaction in the rest frame of the nucleus is ic ~ l/lrriNXBj), where m-Ar is the 
nucleon mass and the Bjorken xbj is the fractional momentum that the struck parton carries. 
At sufficiently low xbj, the distribution is dominated by gluons and the coherence length will 
exceed the nuclear diameter ~ 2A^^^. Modification of parton distributions due to the coherence 
can be studied within the framework of Glauber multiple scattering in the rest frame of a nucleus 
[HH] . More intuitively, saturation phenomena can be studied as multiple parton interactions in 
the infinite-momentum frame. For a hard process with momentum transfer Q, all gluons within 
transverse area l/Q"^ will participate coherently in the interaction. Denoting the nuclear gluon 
structure function as xGa{x,Q'^), the density of gluons in the transverse plane is 



71 K, 



where Ga{x, Q'^) ~ AG{x, Q'^) and G(x, Q^) is the gluon distribution in a nucleon. 

For gluon scattering with cross section cr ~ vra^/Q^, crpA represents the probability of multiple 
gluon scattering. At high or apA <C 1, the system can be considered to be dilute, and the 
perturbative QCD parton model applies. However, for low or apA 3> 1, the target looks 
black to the probe and the saturation regime is reached. The boundary where apA — 1 defines 
the saturation scale Qs- The physical process that leads to saturation is the nonlinear gluon 
fusion gg g, which competes with the gluon emission process g gg. The emission process 
increases the gluon number with rising according to the normal DGLAP evolution, while 
gluon fusion, which is proportional to apA, reduces the gluon number. Saturation occurs when 
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the two processes offset each other and the saturation scale is determined by the self-consistent 
solution of the equation [H3 Ell 

= ^^c.siQs')xG{x, Qs')^, (25) 
- 1 TiR^ 

where the gluon distribution xG{x,Qs'^) is evaluated at a; = 2Qs/y/s. In the saturation regime, 
xGa{x^Qs^) oc l/a{Qs^) and Qs^ oc A^^^. For large enough Qs, the strong coupling constant 
will be small while the density is high, enabling treatment of the non-linear QCD dynamics 
by classical weak coupling methods. This provides the foundation for semi-classical treatment 
of the gluon distribution inside large nuclei jHSlEZI- Often referred to as Colored Glass Con- 
densate (CGC) model, this approach approximates the gluon distribution in large nuclei as 
the Weiszacker- Williams distribution from the classical Yang-Mills field of randomly distributed 
color charges. For a recent review see |17j . 

The initial saturation phenomenon is generic and is independent of the type of hadrons or 
nuclei being collided. However, the growth of Qs^ as A^^'^ suggests that saturation phenomena 
may occur at higher x (or equivalently, lower -y/i) in collisions of heavy nuclei than in p + p 
collisions. Deep inelastic scattering data from HERA indicate that Qs scales as[ 



Qs\x)=Ql(^^) , (26) 

with A ~ 0.2 — 0.3. The saturation scale in heavy ion collisions at RHIC is estimated to be 
Qg^ ^ 2 GeV^ at X ~ 0.02 [HI], roughly the value of Qg^ at x ~ 10""^ in p + p collisions at the 
same energy. 

In contrast to the EKRT final-state saturation model, the initial-state saturation model 
assumes that final gluon production is dominated by gluons below the saturation scale [^2] and 
ignores gluons above the saturation scale, whose yield falls as l/pr^- The density of gluons per 
unit area and unit rapidity produced in the collision is then given by [011 

^ = c Q ' (27) 

Integrating over the transverse area and further assuming the hadronization coefficient is unity 
(one produced gluon becomes one final charged hadron), the observed charged hadron multi- 
plicity in the most central Au + Au collisions is fitted to obtain the gluon liberation coefficient 
c = 1.23 ± 0.20 [HI]. A similar value for c results from numerical calculation of initial parton 
production within the CGC model |100j . 

With these ingredients, the classical weak coupling treatment of initial-state saturation gives 
a prediction for the centrality dependence of the multiplicity density per participant pair [91]: 

-j^dN/df] ^ 0.82 log (^^] , (28) 

where Aqc'd=200 MeV. The centrality dependence results from the variation of Qs^ ~ Ppart with 
the impact parameter according to Eq. (|^ . where 2A/7iR\ is replaced by ppart as the transverse 
density of participant nucleons at fixed impact-parameter pT] . 

The collision energy and rapidity dependence of the multiplicity is governed by Eq. ^I^ . Since 
rapidity ?/~log(l/x), the rapidity dependence of Qs at fixed ^/s is Qs'^{±y) = Qs^{y = 0)6^^^^. 
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Figure 11: y^-dependence of charged particle density per participant pair at midrapidity, for 
central collisions of heavy nuclei {A ~ 200, open and filled points) and p + p collisions (points 
joined by line). Figure from PHOBOS |101j . 



In other words, at forward rapidity one nucleus moves deeper into the saturation regime (larger 
Qs) while the other moves towards the low density domain. The complete expression for the 
multiplicity density in the initial state saturation model is 
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(29) 



where Qs^{s) = Qs'^{so){s/so)'^^'^. This expression contains two free parameters c and Qs'^{so), 
which are fixed at one energy, rapidity and centrality. 



4.2 Multiparticle Production 

Total charged hadron multiplicities were the first published experimental data at RHIC |84j . 
Fig. ^2 shows the energy dependence of the charged particle density at mid-rapidity normalized 
per participant pair, for collisions of heavy nuclei [A ~ 200) and p + p. The nuclear collision 
data are from the AGS, SPS and RHIC (the three highest energy points are from RHIC). The 
energy dependence is smooth, with only logarithmic dependence on ^/s. The growth for nuclear 
collisions is nevertheless faster than that in p + p collisions, qualitatively in agreement with the 
expectation of a larger minijet contribution at higher energy pU]. 

A more differential view of bulk particle production is given by Fig. ^1 which shows the 
centrality dependence of the multiplicity density per participant pair compared to several model 
calculations. In the left panel, the two-component fit and initial state saturation model are seen 
to describe within experimental uncertainties both the centrality dependence and the growth in 
multiplicity with energy. In the right panel, the EKRT prediction apparently fails to describe 
the centrality dependence, while the centrality dependence of HIJING is consistent with the 
data but the normalization is about 10% too low except for p + p collisions. The almost linear 
centrality dependence of the HIJING result is not characteristic of a two-component model. This 
effect may be caused by coherent string fragmentation of the minijet hadronization, which could 
modify the binary scaling of the number of hadrons from jet fragmentation. The two-component 
minijet model jl02j . which assumes independent fragmentation, describes the data well (shaded 
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Figure 12: Centrality dependence of the charged hadron central rapidity density per participant 
nucleon pair for Au + Au collisions from PHOBOS H HUSl and PHENIX ^U^. 200 GeV p + p 
data (leftmost points) from UA5 1051 1106] . Left: theory calculations from two-component mini- 
jet model (shaded bands), two-parameter fit (Eq. (jUJ) (dot-dashed lines) and parton saturation 
model (solid lines) |HS1- Right: data from 130 GeV collisions; theory calculations from Hijing 
HOlini, EKRT saturation (Eq. and two-parameter fit(Eq. 

bands, left panel). Both the HlJlNG model and the two-component minijet model have parton 
shadowing that depends on impact-parameter. Similarly, the saturation models have an impact 
parameter-dependent saturation scale. 

Note, however, that the calculated quantity {Npart) appearing in both the ordinate and 
abscissa in Fig. El is derived using the Monte Carlo Glauber model for the data in both 
panels. The uncertainty inherent in this procedure is demonstrated in Fig. which shows the 
same STAR data in both panels but with (Npart) calculated via the Optical (left) and Monte 
Carlo (right) approach to the Glauber calculation j73j. With the Optical Glauber approach 
the logarithmic growth at low Npart, the hallmark of initial state saturation (Eq. OH)) , is not 
seen, while the final state saturation model (EKRT) still disagrees with the data, though less 
significantly than with the Monte Carlo approach. Though other measurements favor the Monte 
Carlo Glauber (Fig. |HI), the saturation model curves in Fig. [Tulare calculated using Optical 
Glauber, obscuring somewhat the direct comparison to data. The centrality dependence of 
the multiplicity density and its comparison to models therefore requires further clarification, 
but it is apparent that multiparticle production in central collisions can be well described by a 
broad range of theoretical approaches, with about 20% uncertainty both in theory and in the 
measurement of the centrality dependence. All of these models point to initial conditions with 
high initial gluon density. 

4.3 Pseudo-rapidity Distributions 

Fig. 1131 upper panels, show the charged particle pseudorapidity distribution over the full RHIC 
phase space for Au-|-Au collisions at several centralities and energies, measured by PHOBOS 
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Figure 13: dNch/dri/{{Npart)/2) utilizing Optical (left) and MC (right) Glauber calculations 
with the same STAR data [73]. The error bars for adjacent data points are highly correlated. 
Kharzeev-Nardi refers to a two component soft/hard fit [Eq. 



|lU7j . The distributions exhibit two general features: a plateau about midrapidity which broad- 
ens with increasing collision energy, and a forward region whose width is approximately invariant. 
The total charged multiplicity for central collisions at y/s^=200 GeV is 5060 ± 250 |107j . in- 
dicating a qualitatively new regime of accelerator-based experimentation in high energy and 
nuclear physics. 

Figure El lower panels, compare the 200 GeV data from BRAHMS |108j to both saturation 
model jH3] and AMPT [^21 calculations. AMPT combines the HIJING model with partonic and 
hadronic final-state rescattering and includes an alternative baryon pair production mechanism. 
Good agreement with data is found for both model calculations. 

Study of bulk particle production in p{d) +A collisions may help to separate initial from final 
state effects, since a dense medium is not formed in the central rapidity region in such collisions. 
Figure IT31 shows the charged multiplicity density distribution for minimum bias d + Au collisions 
at a/s=200 GeV |109j (see also BRAHMS |110j ). This measurement offers a unique probe of 
particle production, because of the large projectile asymmetry. All panels show the same data, 
which are compared to various model calculations. The data indeed exhibit an asymmetry 
in particle production, with larger multiplicity density towards the direction of the Au beam 
(?7<0). The left panel of the figure compares the data to predictions from the HIJING [701 
and AMPT models The predictions agree with the data except at large negative rapidity, 
the region of the Aw-nucleus fragmentation. Evidently the final state interactions and baryon 
production mechanisms in AMPT provide a better description of this region. 

Using the weak coupling approach to the saturation regime in the initial state saturation 
model, the rapidity dependendence of the charged multiplicity in d + Au collisions is given by 
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where S is the the interaction cross section. Qs,max and Qs,min denote the larger and smaller of 
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Figure 14: Charged hadron distributions over the full RHIC phase space for Au + Au collisions. 
Top: Centrality and y/s^ dependence of dNch/drj from PHOBOS |lU7j . Bottom: Centrality 
dependence of dNch/drj at y/s = 200 GeV, from BRAHMS (circles) |108j . Stars are p + p data 
scaled by Npartf^ |lU5j . Calculations are saturation model (solid) jHS] and AMPT (dashed) 
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Figure 15: dNch/drj for minimum bias d+Au collisions at y^=200 GeV, from PHOBOS 
All panels show the same data. The left panel compares to model calculations from AMPT [221 
and Hijing [701 [^. The middle and right panels compare to saturation model calculations 
(Eq. [20]) using Optical (middle) and Monte Carlo (right) Glauber calculations and slightly 
different proton saturation parameters. 

the saturation scales in the deuteron and the Au nucleus, which vary with rapidity y ~ log(l/x) 
as [Eq. 121] • Qs is assumed to be the same for the deuteron and proton. Since SQs^ ~ Npart 
(H3] . dN/dy ~ Npart{Au) in the Aw-fragmentation region and dN/dy ~ Npart{,d) in the deuteron- 
fragmentation region, replicating the scaling of the phenomenological Wounded Nucleon Model 
p9^. The factor C in Eq. IHOlis determined from the midrapidity charged hadron density in 130 
GeV Au+Au collisions. 

The solid line in the middle panel of Fig. IT^ shows Eq. (j3U|) utilizing an Optical Glauber 
calculation |lllj . (We will not discuss the curve labelled "RQMD".) The calculation overesti- 
mates the measured distribution in the Au-fragmentation region and underestimates it in the 
(i-fragmentation region. A revised version of the calculation (erratum to [.1 1 Ij ) utilizing a Monte 
Carlo rather than Optical Glauber calculation and a slightly different proton saturation param- 
eter achieves good agreement with the measurement (right panel). 

The summary of this section is similar to that of the previous section: both pQCD-based 
models and models incorporating initial state saturation reproduce the pseudo-rapidity distri- 
butions quite well over very broad phase space. The common feature of these models is again 
that the initial energy density is very high, but evidently these observables are not sufficiently 
discriminating to distinguish between the rather different production mechanisms of the models. 
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Figure 16: Rapidity dependence of p^- integrated particle yields from central Au + Au collisions 
at 200 GeV, from BRAHMS. Left: 7r±, K±, p and pjHTl HIHl HIHl- Dashed lines are Gaussian 
fits. Lower panel shows rapidity dependence of {pr)- Right: net protons (difference of p and p 
yields) compared to lower energy collisions jll8j . 



4.4 Rapidity Distributions and Baryon Stopping 

Some time ago, Bjorken postulated that the rapidity distribution at very high collision energy 
should develop a plateau in the central rapidity region, which results from a reaction volume 
that is invariant under longitudinal boost |112j . This assumption leads to considerable simplifi- 
cation of the hydrodynamic equations and is common in theoretical treatments of mid-rapidity 
observables |113j (Sect. 15. Ij) . The plateau in pseudo-rapidity density seen in Fig. HH suggests 
that the fireball near mid-rapidity may indeed be boost-invariant. However, pseudorapidity r] 
only approximates rapidity y. Figure UHl left panel, shows the rapidity dependence of particle 
production separately for pions, kaons and protons. While the rapidity distributions are indeed 
broad, no plateau is observed and except for protons they are Gaussian in shape. Less marked 
but still significant rapidity dependence is also seen for {pt) ■ These distributions may neverthe- 
less result from boost-invariant initial conditions of limited extent in rapidity (e.g. |114t Ill5] ). 
Note that even at LHC energies (ys^=5500 GeV) the initial energy density computed from 
pQCD with saturation scale Psat='^ GeV is not uniform in rapiditv jlltij . so that a boost-invariant 
initial condition may in any case not be the correct high energy limit. 

Baryon number is conserved in the collision and its rapidity distribution should be very 
different from that of produced particles. Since rapidity is logarithmic in energy, it is not 
changed significantly by rescattering. The net baryon rapidity distribution observed in the final 
state is therefore established to a large extent early in the collision, reflecting the mechanisms 
of energy transfer from the colliding nuclei to the fireball. In thermodynamic terms, finite net 
baryon number results in finite baryochemical potential hb (Sect. 15. 4j) . While the distribution of 
total baryon number is difficult to access experimentally, it can be deduced from the net proton 
distribution p — p (Fig. E[ right panel), which is seen to be small but finite at midrapidity. The 
BRAHMS data indicate a net baryon density dN/dy ^ 10 for central Au+Au collisions at y = 
[118] , in agreement with other analyses j761 1120j . 

The mean rapidity loss of leading baryons in lower energy fixed target interactions of protons 



28 



with heavy nuclei is (Ay) ~ 2.5, in contrast to (Ay) ~ 1 for proton-hydrogen interactions |121j . 
The rapidity distribution of net protons in Fig. ^[ taken together with conservation of net 
baryon number, significantly constrains the full net baryon rapidity distribution, resulting in 
{Ay) ~ 2.0 ± 0.2 |118j and a mean energy loss per participating nucleon of AE = 72 ± 6 GeV 
in central Au + Au collisions |118j . In other words, for interacting nucleons about 70% of the 
incoming beam energy of 100 GeV per nucleon is delivered to the fireball. Central collisions 
have about 350 participants, giving ~ 25 TeV transferred from the incoming projectiles to final 
particle production. It is notable both that the net baryon density at midrapidity is small for 
central collisions relative to the 2 x 197 nucleons brought into the collision and that it is finite, 
indicating transfer of baryon number over 5.5 rapidity units. 

The conventional mechanism for baryon transport in hadronic and nuclear collisions is the 
fragmentation of quark-diquark (q — qq) strings |122j . However, such string fragmentation models 
underpredict the baryon stopping measured in nuclear collisions both at SPS |123j and at RHIC 
energies |124j . An alternative scenario considers baryon structure comprising the gauge junction 
which carries the baryon number of three quarks in their fundamental representation |125| I12(jj . 
An implementation of this baryon junction mechanism in HIJING/BB [127 describes well the 
measured baryon stopping at RHIC |124j . Other modified string fragmentation models |128|I129] 
and diquark rescattering [130[ I131j can also provide stronger baryon stopping power than the 
conventional string fragmentation. However, it is not clear at this point which of the baryon 
transport mechanisms are dominant in high energy heavy-ion collisions. 

4.5 Transverse Energy and Energy Density 

Significant transverse energy can only be generated during the collision, through the initial 
interactions of partons from the projectiles and the successive interactions among the produced 
partons and hadrons. Experimentally, E^ is defined as Ej- = J2i EiSin{di), where i sums over 
all final state hadrons. The hadron energy Ei is corrected for conserved baryon number, and 
6i is the angle relative to the beam direction. Due to the dynamics of the expansion, Et in a 
limited rapidity interval will evolve through the lifetime of the collision. In the framework of 
hydrodynamics this dependence is |1321 1133j : 



where tq is the equilibration time. Local thermodynamic equilibrium is characterized by energy 
density e, pressure p, and speed of sound Cq = dp/de. If equilibrium is established at tq and 
maintained throughout the expansion with constant Cg, then 6 = Cq and the observed Et is 
substantially reduced relative to the initially generated Et due to the pAV work performed 
during the expansion |132j . Alternatively, if the system falls out of equilibrium quickly into a 
free-streaming gas, 6 = and there will be no pAV work performed during the evolution, so 
that Et will remain constant throughout the expansion. Final state saturation effects (Sec. 14. 1|) 
reduce early pressure, delaying the onset of hydrodynamic behavior and leading to significant 
reduction in the observed Et |133j . 

It is difficult to disentangle these competing mechanisms based solely consideration of Et 
distributions, but the systematic study of Et together with that of other bulk observables may 
isolate the contribution of longitudinal work. If hydrodynamic flow can be shown to set in early 
from other considerations, the measured Et will provide a lower limit to the initially produced 
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Figure 17: Transverse energy Et for Au + Au collisions at 



=130 GeV, from PHENIX 



compared to Pb + Pb collisions at ys^=17 GeV |135j . Left: Et per participant. Right: E^ per 
charged particle. The Pb + Pb data in both panels have an overall normalization uncertainty of 
±20%, not shown. 



transverse energy at the time of equilibration and thus provide an estimate of the initial energy 
density. 

Fig. El from PHENIX |134j shows the centrality dependence of Et for y^s^=130 GeV 
Au + Au collisions, compared to Pb + Pb collisions at y/s^=17 GeV |135j . In the left panel, Et 
per participant is seen to increase with increasing ^/s for all centralities. As shown in the right 
panel, however, Et per charged particle is largely independent of collision energy, meaning that 
the dependence of Et on energy and centrality closely parallels that of the charged multiplicity in 
Figs. IT^lIT^ Since the average (pt) for charged hadrons in p+p{p) collisions increases significantly 
with ^/s (l^Hj . a constant Et per hadron in heavy- ion collisions implies the existence of pAV 
work due to hydrodynamic expansion. Indeed, hydrodynamic calculations assuming onset of 
equilibration at tqKI fm/c are able to reproduce approximately the centrality dependence of Et 
per charged particle |113j . 

In the Bjorken picture of boost-invariant free-streaming expansion, the initial energy density 
can be expressed in terms of the observed Et |112j : 

- ^^^^ 

The formation time is taken to be tq ~ 1 fm/c and initial system size R ~ 1.2^4^/^, equal to the 
nuclear radius. This expression relies on the assumption that no pAV work is done during the 
expansion |132j and it therefore represents a lower bound to the initial energy density within 
the hydrodynamic framework. PHENIX has measured dET/drj ~ 540 GeV |134] at midrapidity 
for central Au + Au collisions at ys^=130 GeV, resulting in eBj=4.6 GeV/fm^. In comparison, 
NA49 estimates enj ~ 3 GeV/fm^ for central Pb + Pb colhsions at ^/s;^=17.2 GeV^]\. These 
values are provocative: they lie well above the deconfinement energy density predicted by Lattice 
QCD calculations (Sect. 12. 1|) . However, the calculation is based on the assumption rather 
than the demonstration that the onset of hydrodynamic expansion occurs at Tq ~ 1 fm/c. In 
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the following sections we will address this question through measurements sensitive to local 
equilibration at the early, hot and dense phase of the collision, in particular elliptic flow. 

5 Collective Phenomena 

A central question at RHIC is the extent to which the quanta produced in the collision interact 
and thermalize. Nuclear collisions generate enormous multiplicity and transverse energy, but in 
what sense does the collision generate matter in local equilibrium which can be characterized by 
the thermodynamic parameters temperature, pressure, and energy density? Only if thermaliza- 
tion has been established can more detailed questions be asked about the equation of state of 
the matter. 

The initial energy density, whether equilibrated or not, will have strong spatial gradients due 
to the geometry of the colliding nuclei and the dynamics of the collision. Reinteractions among 
the fireball constituents will convert these density gradients into pressure gradients, resulting in 
collective flow of the matter. Collective flow is thus a generic consequence of reinteractions, which 
also lead to thermalization. It is however not sufficient merely to observe collective flow, which 
may be generated both early through partonic reinteractions and later through interactions in 
the dense hadronic gas. It is partonic thermalization and the partonic Equation of State that 
are of interest, but their signals may be masked by the hadronically generated flow that must 
be understood and unraveled. 

In this section we discuss hadronic observables that are sensitive to collective flow and the 
degree of thermalization in nuclear collisions at RHIC. Of particular importance is the azimuthal 
anisotropy of the final hadron spectra in non-central collisions ( "elliptic flow" ) that results from 
the conversion of the initial coordinate-space asymmetry to momentum space via collective 
expansion. The particle mass dependence of flow is an especially sensitive observable, since a 
common velocity distribution for fluid cells radiating particles of different mass will result in a 
characteristic mass dependence of the momentum spectra. 

5.1 Relativist ic Hydrodynamics 

Relativistic hydrodynamics provides the theoretical framework to study collective behavior in 
high energy collisions, with the first such attempts dating back to Landau [138J. We sketch here 
the basic ideas and compare hydrodjTiamic calculations to a wide range of RHIC data. Generally 
good agreement is achieved (though with notable exceptions), providing strong evidence that 
local equilibrium is established early in the evolution of the fireball (r<l fm/c) and that the 
system evolves in accordance with ideal hydrodynamics. Some sensitivity to the equation of 
state is observed, with preference for a deconfined phase early in the evolution. Detailed reviews 
of relativistic hydrodynamics with applications to RHIC data can be found in |113| I139j . 

We first present a simple estimate to assess whether hydrodynamics is a reasonable approach 
to modeling the dynamics of a deconfined phase |139j . Consider a two-fiavor QGP at temperature 
T ~ 200 MeV, which has partonic density n ~ 4 fm~^. Assuming the Debye screening mass 
H = gT as the typical momentum transfer in gluon-gluon scattering, the pQCD cross section 
agg^gg ~ 3 mb glvcs a mean free path A = 1/an ~ 0.8 fm. The time between collisions is 
therefore an order of magnitude smaller than the expected system lifetime of a few fm/c, so that 
the multiple reinteractions necessary for thermalization may occur. 

Hydrodynamic behavior can set in only at a finite time after the collision, when the produced 
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quanta have interacted and relaxed into local equilibrium. Hydrodynamics therefore does not 
address the earliest moments of the fireball evolution, and its initial conditions (density distri- 
butions and flow velocities) must be imposed on the basis of other considerations. The initial 
conditions (entropy, energy, and net baryon number density) are constrained by comparing to 
experimental data such as hadron multiplicities and transverse energy production. 

The hydrodynamic evolution terminates when the system has expanded and cooled to a 
degree that the mean free path exceeds system size and local equilibration can no longer be 
maintained ("freezeout"). Generically, two stages of freezeout are expected: chemical fieezeout 
occurs when the mean free path for inelastic collisions exceeds the system size, whereas kinetic 
freezeout occurs at a later time and a lower temperature, when the elastic mean free path 
also exceeds the system size. The produced hadrons, dominantly soft, are created continuously 
at the dilute periphery of the fireball, according to a specific prescription of kinetic freezeout. 
Though these soft hadrons do not directly transmit signals from the hot and dense early stage 
of the collision, the systematic study of the transverse momentum and mass dependence of soft 
hadron production can provide substantial evidence for early pressure build-up and therefore 
equilibration. Sufficiently detailed and precise comparison of data and calculations will also be 
able to constrain the initial conditions and the EOS at the early stage. 

For ideal, non-dissipative hydrodynamics, the energy- momentum tensor T^'^(x) in the global 
reference frame for a fluid cell at space-time coordinate x is given by |113j 

T^'^{x) = [e{x) + p{x)]u^{x)u''{x) - p{x)g^'\ (33) 

where e(x) is the energy density, p{x) is the pressure, and M^(a;) is the four-velocity of the cell. 
Correction for non-ideal hydrodynamics adds a term that is the product of the shear viscosity rj 
with the thermally averaged gradient of the velocity field |14Uj . 

The equations of motion result from local conservation of energy and momentum, 

9^T^'^(x) = 0(z/ = 0,...,3). (34) 

Additional equations result from the conservation of M different charges (net baryon number, 
net strangeness, electric charge), 

d,j^{x)=0{z = l...,M), (35) 

where jf (x) = ni{x)u^{x) is the current density in the global frame and ni{x) is the local charge 
density. 

Expressions (jH^ and (jH^j) comprise 4 + M differential equations for 5 + M fields: the three 
components of the flow velocity, the energy density, the pressure, and the M charge densities. 
The system of equations is closed by the equation of state (EOS) p{e,ni), which relates the 
pressure, energy density, and conserved charge densities. Most applications of hydrodynamics 
to RHIC data use a similar structure for the equation of state |139j : a plasma phase of massless 
partons with a bag constant, a hadronic phase consisting of a gas of free hadrons and resonances, 
and a flrst order phase transition with a large latent heat connecting the two phases. 

The initial conditions at the onset of hydrodynamic expansion must be specifled from external 
input, usually either the entropy or energy density, with distribution in the transverse plane 
according to that for binary collisions or participants nucleons |139j . Saturation initial conditions 
have also been considered |141j . The conventional implementation of freezeout is via the Cooper- 
Frye prescription |142j which corresponds to an instantaneous transition from zero to inflnite 
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Figure 18: Radial dependence of the flow velocity from full hydrodynamic calculations of cen- 
tral Au + Au collisions at RHIC. Left: at various proper times Right: dependence at 
hadronization on equation of state |145j . The left panel is most comparable to LH8. 



mean free path, i.e. from ideal hydrodynamics to free streaming. The spectrum of hadron 
species i at freezeout is given by 

^dNi Qi f 1 



E—^ = / — p'^dV., (36) 

where the integral is carried out over the hypersurface on which the freezeout conditions 
are met. ^i{x) is the local chemical potential for species i and T(x) is the local temperature. 

Having specified the EOS and initial conditions, the differential equations ()34|) and ()35|) 
are integrated numerically to freezeout, where the stable hadrons and resonances are generated 
according to Eq. ()36|). Integration of the full three dimensional hydrodynamic equations is 
a daunting task |143j . The assumption of longitudinal boost invariance is often made, with 
the imposed symmetry reducing the number of coupled equations and simplifying the numerical 
problem considerablv |113j . This approach is however only applicable to mid-rapidity observables. 

The instantaneous freezeout embodied in the Cooper-Frye prescription is unphysical. A more 
realistic though calculationally more intensive transition to on-shell hadrons results from coupling 
the hydrodynamic evolution to a kinetic transport model |144t I145j . However, at present the 
experimentally accessible observables exhibit no significant variation between this approach and 
the simpler Cooper-Frye algorithm |113j . 

Fig. ITHl shows the radial profile of the transverse flow velocity Vr resulting from two different 
hydrodynamic calculations for central Au + Au collisions at RHIC |113l 1145] . The left panel 
shows the time dependence of Vr with an EOS incorporating a first order phase transition. In 
the bulk (r< ~ 6 fm), the buildup of transverse velocity due to pressure is rapid, achieving a 
roughly linear gradient ~ 0.07//m that persists for the lifetime of the fireball. The velocity near 
the dilute surface varies strongly with radius at early times due to the initialization of matter 
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Figure 19: Identified particle inclusive spectra for centrality-selected Au + Au collisions at 
.^s;;;;=200 GeV. Left: 7r±, K±, p and p from PHENIX j^. Right: A + A and K'^, from STAR 
[TIHllTTfj . 



in the mixed phase, which has vanishing pressure gradient, and in the hadronic phase at the 
largest radii. This surface feature is eventually overtaken by the expanding plasma at higher 
pressure. The largest system size is achieved at r ~ 10 fm/c; by 15 fm/c the freezeout surface 
is contracting inwards. 

The right panel of Fig. shows the transverse rapidity yx = tanh^^ Vr for three different 
equations of state |145j : a hadronic resonance gas (RG), plasma and hadronic phases linked by 
a mixed phase with latent heat 0.8 GeV/fm^ (LH8), and mixed and hadronic phases only (i.e. 
infinite latent heat, LHoo). The flow profiles are shown for constant energy density e = 0.45 
GeV/fm^, where in this calculation the hydrodynamic evolution is terminated and the produced 
hadrons propagated further using a kinetic transport model. The resonance gas EOS is seen to 
be quite stiff, generating a higher transverse velocity gradient than those containing a mixed 
phase. The presence of the plasma phase also provides significant pressure, generating twice the 
flow velocity at large radius than the case where it is absent (LH8 vs LHoo). Similar to the 
calculations in the left panel, LH8 also produces a radial velocity gradient ~ 0.07/fm late in 
the evolution. Note that constant energy density does not correspond to constant proper time 
in this calculation. The double-valued loop at large radius for LH8 is also due to matter on 
the dilute surface initially generated in the mixed or hadronic phase, freezing out rapidly after 
modest radial expansion. 



5.2 Transverse Radial Flow 

In the hydrodynamic picture just described, the observed final state hadrons freeze out from 
fluid cells that are in local equilibrium but that have finite transverse velocity relative to the 
lab frame. Since the thermal sources for all hadron species are boosted with the same velocity 
distribution, the hydrodynamic expansion should result in a characteristic mass dependence of 
the transverse momentum spectra for momenta on the order of the particle mass. 

We first look at the systematic features of the data. Fig. shows inclusive transverse 
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momentum spectra for tt^, K"*^, p and p from PHENIX [TH] (left panel) and A + A and from 
STAR jl46l I147j (right panel), for centrality-selected Au + Au collisions at ys^=200 GeV. 
The shape of the baryon spectrum changes qualitatively from peripheral to central collisions. 
Relative to the meson yields, the baryon yields in central collisions are suppressed at low pt 
and enhanced at higher p^, with a marked change of slope at 1-2 GeV/c. For pt>'2 GeV/c, the 
baryon yields exceed the meson yields (Sect. 16. 5p . 

Fig. 1201 compares hydrodynamic calculations |148tlll3] to measured pt spectra from Au + Au 
collisions of vr", K+ and p |149|ll2Utll5Uj and fl~ jl51j . For the upper left, upper right, and lower 
left panels, the parameters of the calculation were fixed by fitting the tt"*" and p distributions 
in central collisions, resulting in an equilibration time r = 0.6 fm/c with temperature T = 340 
MeV and energy density e = 25 GeV/fm^ at the core of the fireball. The remaining curves 
in those panels are then predictions of the model. Overall agreement with the data is good. 
Deviations are seen at low p^ for the pions, now understood to be due to the imposition of 
chemical equilibrium through to kinetic freezeout |113j . More significant disagreements are seen 
for Pt>2 GeV/c in the most peripheral collisions. This may delineate the region of applicability 
of the hydrodynamic approach, since the fireball is smallest for peripheral collisions and high p^ 
particles require the greatest number of collisions to thermalize jll3j . 

The lower right panel of Fig. 1201 compares the pt spectrum of Q~ in central 200 GeV Au + Au 
collisions to a similar calculation, adjusted for the higher collision energy and with a chemical 
non-equilibration EOS in the hadronic phase |152j . The steeper set of curves results from 
decoupling at energy density e = 0.45 GeV/fm^, corresponding to kinetic freezeout at the onset 
of the hadronic phase. The shallower set of curves results from decoupling at the same time as 
the pions at e = 0.075 GeV/fm'^, so that the Q~ receives the full boost from the hadronic phase. 
(Solid lines are for zero initial flow, dashed lines are for small but finite radial flow at the onset 
of hydrodynamic behavior |ll52j ). It has been proposed that multistrange baryons decouple from 
the flow much earlier than non-strange hadrons due to absence of strong resonances with pions 
|153j . Evidently, purely partonic hydrodynamic flow (lower curves) does not generate sufficient 
transverse velocity to describe the Q~ spectrum and there is significant contribution from the 
hadronic phase |113j . 

The systematic behavior of fiow-related observables is extracted by fitting measured spectra 
to the phenomenological Blast Wave parameterization |154| 1155] , which results from modeling 
the system at freezeout as an ensemble of transversely boosted Boltzmann distributions. The 
transverse velocity distribution is a parameterization of the radial dependence of the fiuid cell 
velocity distribution at freezeout from the full hydrodynamic calculation (Fig. ITH|) : 

Mr) = /3.(-^)", (37) 

where R is the radius at freezeout and f3s is the transverse fiow velocity at the surface, n = 1 
reasonably approximates of the full calculation. Assuming that kinetic freezeout occurs instan- 
taneously at all radii, the hadronic spectra are given by |155j : 

dN fR (pTsinhp\ /m^coshpN 

(X rdr mrlo Ki , (38) 

mTdrriT Jo \ I J \ I J 

where ttit = \Jpt^ + vd^ ^ pij) = tanh~^/?T(^) is the transverse rapidity and T is the local 
temperature. Non-instantaneous freezeout results in additional terms in the integrand which 
can change the spectrum shape significantly |155j . 
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Figure 20: Transverse momentum spectra of identified hadrons from Au + Au collisions at 
RHIC, compared to liydrodynamic calculations. Upper left, upper right, and lower left: 130 
GeV Au + Au data from PHENIX [Hn] and STAR 11201 IM, calculations from [UHl. Lower 
right: 200 GeV central Au + Au data from STAR [l5T], calculations from [113J. 
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Figure 21: Freezeout temperature Tfo and mean flow velocity {/St) for centrality-selected 130 
GeV Au + Au collisions, from PHENIX [156|. Filled points are Blast Wave flts to data, open 
boxes are from a full hydrodynamic calculation. The dotted line is for blast wave with centrality- 
independent T/o=128 MeV. [THTj . 



Figure |^ shows the collision centrality dependence of the freezeout temperature Tfo and 
mean flow velocity {Pt) for a Blast Wave flt to the measured vr, K and p spectra for 130 GeV 
Au + Au collisions |156j . and for a full hydrodynamic calculation similar to that in Fig |20l P-57j . 
The dashed line shows (Pt) for flxed Tfo=l2S MeV, the value used in the full calculation, allowing 
a direct comparison of the parameterization with the theory. The agreement between the Blast 
Wave parameterization and the full calculation is good for more central collisions, indicating 
that the Blast Wave formulation contains the essential freezeout features of the hydrodynamic 
calculation. All approaches show lower Tfo and higher (Pt) for more central collisions, indicating 
a longer expansion time until freezeout. Marked deviations are seen only for the most peripheral 
collisions, where the region of local thermalization, if any, may be small and short-lived. Figs. 
EUl and 1^ suggest that hydrodynamic flow dominates the flreball evolution for all but the most 
peripheral collisions. 



5.3 Anisotropic Flow 

In non-central nuclear collisions the overlap region is azimuthally anisotropic. Fig. EH left panel, 
shows the density of binary collisions in the transverse plane, a common basis for calculating the 
initial energy density, for Au + Au interactions with impact parameter b=7 fm |113j . For hydro- 
dynamic evolution, the initial spatial anisotropy will generate azimuthally anisotropic pressure 
gradients that are stronger in the reaction plane (horizontal in the flgure) than perpendicular 
to it. This will result in a spatially anisotropic momentum distribution which is experimentally 
observable [158j. 

At mid-rapidity the odd harmonics vanish by symmetry and the leading anisotropy is ellip- 



37 



-10 



0.1 



RHIC 
EOS I 



0.05 



-10 



5 10 

X (fm) 







T (fm/c) 



Figure 22: Azimuthal spatial and momentum anisotropy for non-central (b=7 fm) Au + Au 
collisions. Left: density of binary collisions in the transverse plane |113j . Right: time evolution 
of spatial eccentricity and momentum anisotropy ep J159'. Solid lines are for EOS with phase 
transition, dashed line for massless ideal gas at very high temperature. 



tical. The initial spatial eccentricity is defined as 



(39) 



where (...) indicates the average weighted, for instance, by energy density. The corresponding 
momentum space anisotropy in the hydrodynamic framework is 

^ Jdxdy{T-- + Tyy)- ^ ' 

This momentum space anisotropy ( "elliptic flow" ) results from interactions within the medium 
and therefore develops over time as the fireball evolves. Figure |22l right panel, shows the 
time evolution of both and for b=7 fm Au + Au collisions. The solid lines are from a 
hydrodynamic calculation with a first order phase transition, similar to that in Fig. 1201 |159j . 
The spatial eccentricity ex is large at the onset of hydrodynamic flow (~ 0.27) but decreases 
continuously with time. The pressure is sufficient to drive negative prior to freezeout. In 
contrast, grows rapidly but saturates at r ~ 6 fm/c due to the low pressure in the mixed 
phase. A small increase in is generated in the hadronic phase, but almost all of the final 
momentum asymmetry is generated in the partonic phase at early time. The early buildup of 
momentum anisotropy is seen not only in hydrodynamic calculations but also in kinetic transport 
models |16mil6l] . Elliptic fiow is thus a key observable of collective hydrodynamic behavior and 
thereby thermalization at RHIC. It is predominantly generated early in the evolution and is 
potentially sensitive to the properties of the partonic stage of the collision. 

Experimentally, fiow is measured by fitting the triple differential invariant momentum dis- 
tribution in a Fourier series in azimuthal angle (158^ (see also ^1^ ): 

d^N 1 d'^N / \ 
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Figure 23: Elliptic flow of charged hadrons from 130 GeV Au + Au collisions, from STAR 
[77] . Left: centrality dependence of pT-integrated V2. Central events correspond to rich/nmax ~ 
1. Boxes indicate range of expected values from hydrodynamic calculations jTTj. Right: pt 
dependence of V2 for minimum bias collisions, calculations from |157j . 



where \E'r represents the orientation of the reaction plane in the event. The Fourier coefficient f„ 
measures the asymmetry of order n. By symmetry, f i = at midrapidity and the leading term 
is elliptic flow t>2. \E'r is of course not directly measurable, but it can be estimated based on the 
azimuthal distribution of all measured particles in the event. However, since the multiplicity is 
finite, the resulting flow coefficients Vn must be corrected for the finite resolution of the estimated 
reaction plane orientation |[158j. 

Fig. 1221 compares V2 for charged hadrons from 130 GeV Au + Au collisions ^77] with hy- 
drodynamic calculations. The left panel shows the centrality dependence of f 2 integrated over 
transverse momentum. Good agreement of the expectations from hydrodynamics with the data 
is seen for nch/nmax>0-5, corresponding to 6< ~ 7 fm. The right panel shows the px depen- 
dence of V2 for minimum bias collisions. Since multiplicity is largest in central collisions and 
the anisotropy is largest in peripheral collisions, the greatest contribution to this measurement 
comes from intermediate impact parameters which have both large V2 and signiflcant multiplic- 
ity. The distribution is compared to hydrodynamic calculations with various equations of state 
|157j . The agreement of calculations with data is impressive, though it is notable that expecta- 
tions for unidentifled charged hadron flow are very similar for an EOS with a flrst order phase 
transition (EOS Q, Tjo =120 and 140 MeV) and a purely hadronic resonance gas (EOS H). 
The only calculation that is excluded is the low-density limit (LDL) from a non-hydrodynamic 
kinetic transport approach |157j . 

Fig. shows V2{pt) separately for various identifled particles from STAR |163j . Similar 
to the azimuthally averaged px distributions in Figs. ^1 and 1201 V2{pt) also exhibits a mass 
dependence. Below 2 GeV/c, proton V2 is signiflcantly smaller than pion V2 for a given value of 
Pt- In hydrodynamic terms, the origin of this effect is the same as the flattening of the inclusive 
Pt spectrum for higher mass: the velocity boost depletes the low pt region in favor of higher pt 
[164] . This systematic trend is seen explicitly in the right panel, where f 2 for different identifled 
hadrons splits according to the mass, consistent with hydrodynamic predictions |164j . At higher 
Pti the hadronic species dependence of f 2 changes due to breakdown of the hydrodynamic model 
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Figure 24: V2{pt) for identified hadrons from STAR |163l I146j compared to liydrodynamic 
calculations from [164J. Left: 71"^ and p+p for 130 GeV minimum bias Au + Au, figure from 
[139j . Right: vr, K, p and A from 200 Gev minimum bias Au + Au as compiled in |147j . 



at high pt, a point to which we will return later when discussing hard probes. 

PHOBOS has reported p^-integrated elliptic flow of charged particles over very broad phase 
space (|r7|<5) (165t 1166] . showing a rapid decrease in integrated V2 away from mid-rapidity. A 
fully three-dimensional hydrodynamic calculation |167j agrees with the measurements at mid- 
rapidity but disagrees significantly at high rj, predicting only weak variation of V2 out to r/ ~ 4. 
Hydrodynamic behaviour is thus limited to low pt hadrons (< ~ 2 GeV/c, depending on particle 
species) at mid-rapidity for more central collisions. Outside of these limits the system size and 
lifetime may be too small for full thermalization to develop fTT^. In contrast, a microscopic 
calculation based on a string model incorporating string excitations and hadonic rescattering 
can broadly describe the full phase space distribution of p^-integrated elliptic flow jl68j , though 
detailed comparison of such an approach to pT^-differential flow and its mass dependence has not 
yet been carried out. 

The left panels of Fig. I^also show hydrodynamic calculations for two EOS (H is a hadronic 
resonance gas, Q contains a first order phase transition) and freezeout temperatures 120 and 140 
MeV. Pion V2 exhibits little sensitivity to the EOS or freezeout temperature, whereas proton 
V2 favors EOS Q(120), with a phase transition and the longest evolution. In |145t IllJ] it was 
concluded that a hadronic resonance gas cannot reproduce the mass dependence of elliptic flow, 
which therefore requires a partonic phase. While the difference between the plasma and purely 
hadronic scenarios in Fig. |^is modest and clear discrimination cannot be made based on these 
data, it is nevertheless evident that elliptic flow does have sensitivity to the EOS in the early 
stage of the collision. Detailed model comparisons to the higher precision data in the right panel 
and future data, especially multi-strange baryon and D- meson V2, will sharpen the arguments 
considerably and may provide significant constraints on the equation of state. 

5.4 Statistical Distribution of Hadron Yields 

The previous sections discussed the effect of early equilibration on expansion dynamics and its 
reflection in the systematic behavior of transverse momentum spectra and their anisotropics. 
Equilibration of the fireball may also be evident in the pr-integrated hadronic yields, which 
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Figure 25: Hadronic yield ratios at midrapidity from central Au + Au collisions at 130 and 200 
GeV. Horizontal bars show the statistical distributions resulting from fits to the data |172t I173j . 

should be statistically distributed according to the thermodynamic conditions at chemical freeze- 
out |169t I170l I17H [T7^ I173j . Measured yields at midrapidity from nuclear collisions at RHIC 
are available for a wide variety of hadron species. In this section we compare relative particle 
abundances to expectations from a statistical model. 

The statistical description of a system with many degrees of freedom is formulated using 
the grand canonical (GC) ensemble, in which conservation laws are enforced on the average via 
chemical potentials /z. The GC partition function for a hadron resonance gas at temperature T 
in volume V is jl73 !: 

logZ(T, V,fl)=Y: log^.(r, V, II), (42) 

i 

where i sums over all hadrons with masses less than ~2 GeV/c^ and fl = {fis, f^s^ f^g) are the 
chemical potentials for baryon number, strangeness and charge. For the hadron species carrying 
baryon number Bi, strangeness Si and charge Qi, 

Z,iT, V,fl) = ^ ±p^dp\og{l ± Ke-^^^). (43) 

Here + is for fermions and — is for bosons, is the spin-isospin degeneracy, (3 = 1/T, = 
y^pM-m?, and 

Ai(T, fl) = exp ^ j . (44) 

Repulsive interactions between hadrons and the effects of resonance decay are taken into account [l 71 j . 
The imposition of local strangeness and charge neutrality means that the resulting distributions 
depend only on the temperature T and baryochemical potential /i^. 

Fig. 1231 shows the ratios of a wide variety of hadronic yields measured at midrapidity for 
central Au + Au collisions. The ratios span three orders of magnitude (note the scaling of 
Q/7T~). Also shown are expectations for a statistically distributed population emitted by an 
equilibrated medium with T = 176 MeV, /iB=41 MeV (^/s;;;;=130 GeV, left) or T = 177 
MeV, /iB=29 MeV (200 GeV, right). Agreement between data and model is good: the hadron 
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population is statistically distributed, consistent with the model sketched in Eqs. (jl^ - Pij) . The 
fitted baryochemical potential /i_B is small, indicating low net baryon density in the medium. 
The chemical freeze-out temperature T appears to be limited (i.e. varies little with y/s) at a 
value close to the QCD phase transition temperature from lattice calculations (Sec. 12. These 
phenomena can be accomodated in a picture in which chemical freeze-out occurs in central 
nuclear collisions at the hadronization boundary of a thermalized, deconfined plasma phase, with 
the subsequent evolution of the hadronic gas being moderated by (quasi-)elastic collisions |178j . 
Additional support for this picture is supplied by the systematics of strangeness production, 
which may indicate that strangeness percolates over a much larger volume in nuclear collisions, 
as described by a grand canonical ensemble, than in elementary nucleon-nucleon collisions, which 
require a canonical ensemble description with explicit local strangeness conservation |173j . 

Does the observation of a statistically distributed population of final state hadrons require 
a chemically equilibrated source? A statistically distributed population may result simply from 
phase space dominance |174j . Surprisingly, hadron populations are found to be statistically 
distributed in p + p and even e"*" + e~ collisions |175t I176j . For a final state comprising many 
particles, a multiplicity measurement corresponds to integration over a large phase space volume 
and details of the matrix element for producing any specific state are unimportant. If the 
phase space-averaged matrix elements obey general scaling rules and do not exhibit strong 
correlations or dependence on -^i, then the resulting hadron distributions will populate phase 
space statistically jl74j . Discrimination between phase space dominance, in which T and fiB are 
simply Lagrange multipliers parameterizing the phase space, and emission from a thermal system, 
having physical parameters temperature T and chemical potential fiB, may be achievable via 
multiparticle correlation measurements |174|I177] . Such measurements are only now maturing at 
RHIC (see |178j for further discussions). The current experimental data are however consistent 
with a picture of hadronic freeze-out from a chemically equilibrated source. 

5.5 Space-time Evolution 

We have so far concentrated on momentum-space observables to infer the initial conditions, 
degree of thermalization, and dynamics of the expansion. Direct measurements of the space-time 
evolution of the fireball provide complementary probes of the system dynamics. For instance, 
measurement of a long system lifetime may indicate the presence of a low pressure phase which 
has stalled the expansion, as expected from a first order phase transition with a soft equation of 
state in the mixed phase |179j . 

It was recognized long ago that intensity interferometry of pairs of identical particles is 
sensitive to the geometry of the source, both in astrophysical systems and in elementary particle 
collisions (Hanbury Brown-Twiss or "HBT" interferometry I18X] ). For bosons, the wave- 
function symmetry results in an enhanced coincidence rate for pairs having small momentum 
difference, with the momentum range of the enhancement varying inversely with the space-time 
dimensions of the source. The application of intensity interferometry to high energy nuclear 
collisions provides a unique probe of the dynamic properties of the fireball |182j . The correlation 
function of identical pions encodes the system geometry and expansion dynamics at kinetic 
freezeout, late in the fireball evolution. Its projections relative to the beam direction and to the 
pair mean momentum refiect various aspects of the expansion dynamics and space-time extent of 
the source. Fortunately, the large multiplicities generated in nuclear collisions provide the large 
pair statistics necessary for detailed investigation of the multi-dimensional correlation function 
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The pair correlation function C(q, K) is a function of the pair's relative 4-momentum q = 
Pi —P2 and mean momentum K = +P2), and is related to the Wigner density of the emitting 
source S{x,K)^^: 

= Md^J [mWJ^^ \Id^xS{x,K)f ■ ^^^^ 

S{x, K) can be understood as the probability that the source emits a particle with momentum K 
from space-time point x. Experimentally, C (q, K) is constructed from the ratio of the measured 
pair distribution to a distribution of mixed pairs drawn from different events |184j . The measured 
correlation function C (q, K) has structure due to Bose-Einstein statistics, resulting in enhanced 
probability for small |q|, and to final state interactions, which mask the BE enhancement and 
which must be disentangled to extract geometric quantities. 

Within a Gaussian approximation to the spatial distribution of S{x, K) |182j . the correlation 
function for central collisions is characterized by its projections onto the orthogonal longitudinal, 
outward and sideward directions, which are parallel to the beam, parallel to Kf-, and perpen- 
dicular to Kt respectively. Kt is the pair momentum vector perpendicular to the beam. The 
conjugate radius parameters are Ri, Ro and Rs |185| I18(i| I187j : 

C(q, K) ^ 1 + A exp {-R^qf - Rlql - Rlqf) . (46) 

In high energy nuclear collisions the source expands longitudinally and transversely. In a 
hydrodynamic picture, a pair of identical pions with small momentum difference are unlikely 
to be emitted from different fluid elements with significantly differing velocity. The parameters 
R therefore do not reflect the dimensions of the entire source, but rather the rms widths of 
the effective source that emits particles with momentum Kt ( "regions of homogeneity" ) ^18R] . 
The systematic dependence of R on pair momentum is a rich source of information about the 
expansion dynamics of the fireball |182j . For an infinite Bjorken (boost-invariant) source at 
temperature T the longitudinal radius is 



R',{Kr) ^ rl^, (47) 
though corrections for realistic sources are significant |19nj . An extended lifetime of the system 



may be observable via the combination 

Rl-Rl^PW). (48) 

which is sensitive to the duration of particle emission. Here, = Kt/ Kq is the transverse veloc- 
ity of the pair and (t^) = (t^) — (t)^ is the variance of the particle emission time. Hydrodynamic 
calculations |179j indicate that a large value of R^ — R^, or specifically the ratio i?o/i?s ^ 1, is a 
rather generic indication of a very soft equation of state stalling the expansion, which can only 
arise from the presence of a mixed phase. 

Fig. 1211 shows the measured HBT parameters for central Au + Au collisions at 130 GeV. The 
parameters are not large (<8 fm) relative to the radius of a Aw-nucleus and exhibit negligible 
change from measurements with heavy nuclei at much lower ^/s |191j . Most significantly, the 
ratio Ro/Rs < 1 in 0.2<Kt<1.2 GeV/c [Ml 11121 CHS!, in contrast to the expectation that 
Ro/ Rs ~ 1.5 for a long-lived source. 

The reduction in R with increasing pair Kt is qualitatively consistent with expectations 
from a longitudinally and transversely expanding source. However, quantitative comparison of 
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Figure 26: Charged pion HBT parameters for 130 GeV central Au + Au collisions from STAR 
jl91j and PHENIX |192j . compared to hydrodynamic calculations |113j . Curves are described 
in text. 



the measured radius parameters to the boost-invariant hydrodynamic calculations in the figure 
reveals significant disagreements |113j : neither the magnitude nor the Kt dependence are well 
described. Modifications to the calculation such as earlier freezeout (dotted line), faster buildup 
of flow at the partonic stage (short dashed line), or earlier onset of hydrodynamic behavior 
(long dashed line) either do not fully rectify the disagreements or worsen the agreement with 
inclusive spectra. Likewise, relieving the boost-invariance condition and imposing a more realistic 
treatment of freezeout than the instantaneous Cooper-Frye algorithm do not fully resolve the 
problems. Introduction of finite viscosity generates Ro/Rs ~ 1, but at the expense of significant 
disagreement with V2 measurements |14Uj . A hybrid parton/hadron cascade calculation generates 
Ro/Rs ~ 1 |194j . However, the freezeout in this calculation occurs earliest at small radii, in 
contrast to hydrodynamic calculations where the freezeout surface generically propagates inward 

HE]. 

For non-central collisions, the azimuthal modulation of HBT parameters relative to the re- 
action plane orientation provides a more detailed view of the geometry and dynamics of the 
source at freezeout |195t 119b] . Recent measurements of the azimuthal dependence of HBT radii 
|197j indicate that the source at freezeout is asymmetric and extends out of the reaction plane 
(ea;>0), consistent with a rapid pressure buildup and early freezeout. This picture should be 
contrasted with the calculations shown in Fig. |22l right panel, in which hydrodynamic expansion 
prior to freezeout lasts long enough to turn the spatial anisotropy slightly negative {ex<0, source 
extended in the reaction plane). 

Reconcilation of HBT radius measurements and hydrodynamic calculations has not yet been 
achieved within the currently available theoretical framework (the "HBT puzzle") |113j . How- 
ever, insofar as the HBT correlations are most sensitive to the system properties at kinetic 
freezeout, this may indicate a lack of understanding of the late expansion stage and ultimate 
breakup of the system rather than the dynamics at the earliest, hot and dense phase. The effect 
of the source opacity on the correlation function also remains an open question. 
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6 Hard Probes 



The lifetime of the hot and dense matter produced in heavy ion coUisions at RHIC is estimated 
to be on the order of a few fm/c. Its initial transverse radius is about 6 fm and it undergoes 
rapid longitudinal and transverse expansion. Due to the transient nature of the matter, external 
probes cannot be used to study its properties. Fortunately, the dynamical processes that pro- 
duce the bulk medium also produce energetic particles through hard processes. These energetic 
particles penetrate the bulk matter and reach the detectors as distinct signals. Study of these 
energetic particles and their interaction with the medium, analogous to the method of computed 
tomography (CT) in medical science, will yield critical information about the properties of the 
matter that is impossible to obtain from the soft hadrons produced by hadronization of the bulk 
medium. 

Properties of a medium are conventionally studied using the scattering of particle beams. 
In deeply inelastic scattering (DIS) experiments, for example, leptons scatter off a nucleon 
via photon exchange with quarks. The response or correlation function of the electromagnetic 
currents, 

= l^j d'xe^'^-^A I jr(0)iT(a^) I A) , (49) 

is a direct measurement of the quark distributions in a nucleon or nucleus, where 
j^"^{x) = J2q^q4'qix)'lfj,4'q{x) IS the hadroulc electromagnetic current. Such experiments have 
provided unique information about the partonic structure of nucleons and nuclei and their QCD 
evolution UMllini- 

The scattering technique is not applicable to the dynamic systems produced in heavy-ion 
collisions. However, it has been shown that the thermal average of the above correlation function 
gives the photon emission rate from the evolving system |2UUj . The emission rate depends mainly 
on local temperature or parton density, while the total yield depends on the entire history of the 
system evolution. The properties and dynamics of a strongly interacting system may therefore be 
probed via the measurements of photon and dilepton emission. Additional information is encoded 
in the resonance properties and medium modification of the emitted virtual photons. Screening 
in a color-deconfined medium leads to dissociation of bound states, resulting in quarkonium 
suppression j22]- The color screening arises from the strong interaction between quarks and 
gluons at high density and temperature, which also causes the attenuation of energetic partons 
propagating through the medium. Such an effect underlies the phenomenon of jet quenching and 
the application of jet tomography to probe the dense matter generated in high-energy heavy-ion 
collisions [MIM- 

Charmonium suppression significantly in excess of normal nuclear absorption has indeed 
been observed in fixed-target heavy- ion collisions at the CERN SPS by NA50 [15j, leading to 
speculation that dense, color deconfined matter has been created in central Pb + Pb collisions at 
SPS energy. A direct photon excess over hadronic sources has been observed at the SPS |2U3j . 
and low mass dilepton spectra show signs of medium modifications of hadron properties [TB] . 
Reviews of these results can be found in [TTl I2(J4[ l2U5j . At present, however, the experimental 
investigation of real and virtual photon production at RHIC is just beginning. In this section 
we concentrate rather on the theory and phenomenology of jet quenching at RHIC, for which a 
large body of mature data is available. 

Theoretical studies of medium-induced partonic energy loss date back to an unpublished 
paper by Bjorken, who calculated the energy loss due to elastic scattering in a hot medium. A 
simple estimate is given by the thermally averaged energy transfer Uei ~ Q^/'^'UJ of the jet parton 
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Figure 27: Binary collision-scaled ratio of charged liadron and vr'' inclusive spectra from 200 
GeV Au + Au and d + Au relative to that from p + p collisions, from BRAHMS 208 (upper left), 
PHENIX|2nn] (upper right), PHOBOSgini (lower left) and STARt74j (lower right). See Sect. 
16.41 for details. 



to a thermal parton with energy cu, where q± is the transverse momentum transfer of the elastic 
scattering. The resulting elastic energy loss 



is sensitive to the temperature of the medium but is in general small relative to the radiative 
energy loss discussed below. Here, fi is the Debye screening mass and C2 is the Casimir factor 
of the propagating parton in its fundamental representation. Elastic energy loss can also be 
calculated within finite temperature QCD field theory |207j . with similar results. 

Radiative partonic energy loss was first estimated using the uncertainty principle |211j . The 
first theoretical study of QCD radiative partonic energy loss, by Gyulassy and Wang |212[ 
I213j . modeled multiple parton scattering using a screened Coulomb potential and found that 
Landau- Pomeranchuk-Migdal (LPM) interference effects |214[ I215j play a crucial role. Baier 
et al. (BDMPS) |216j later considered gluon rescattering, which was found to dominate the 
gluon radiation induced by multiple scattering in a dense medium. These initial studies have 
been followed by many more recent works on the subject, including a path integral formulation 
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Figure 28: Left: Dihadron azimuthal correlations at high pt for p + p, central d + Au and central 
Au + Au collisions (background subtracted ) from STAR |223^ l74j. Right: azimuthal correlation 
strength with the reaction plane {v2, Eq. ()4Hl ) of high pt hadrons in non-central Au + Au 
collisions from STAR |224j . Different symbols correspond to different methods for calculating 
V2 (reaction-plane, 2- and 4-particle cumulant; see Sect. I6.7p . 



[217j and an opacity expansion framework |218| I219| I22()j which is suitable for studying multiple 
parton scattering in a thin plasma. The unique feature of radiative energy loss in QCD is its 
non-linear dependence on distance, arising from the non-Abelian LPM interference effect in a 
QCD medium. In this review, we will take the approach of twist-expansion |22H I222j , since 
it connects naturally the discussions of partonic energy loss in a cold nuclear medium and hot 
quark gluon plasma. 

Before continuing with the theoretical discussion of partonic energy loss and its effects in 
heavy ion collisions at RHIC, it is worthwhile first to gain an impression of the reach of the 
available data addressing this physics and the magnitude of the effects under discussion. Figures 
l27l and l28l show the most significant high pt measurements made at RHIC thus far. Both figures 
incorporate measurements of ^/s=200 GeV p+p, d+Au and centrality-selected Au+Au collisions 
at RHIC, with the simpler p + p and d + Au systems providing benchmarks for phenomena seen 
in the more complex Au + Au collisions. 

Figure EH shows the ratio of inclusive hadron yields in Au + Au and d + Au to p + p, scaled 
by (Nhin) to account for trivial geometric effects. A striking phenomenon is seen: large pt 
hadrons in central Au+Au collisions are suppressed by a factor 5 relative to naive expectations. 
Conventional nuclear effects, such as nuclear shadowing of the parton distribution functions and 
initial state multiple scattering, cannot account qualitatively for the suppression. Furthermore, 
the suppression is not seen in d + Au but is unique to Au + Au collisions, proving experimentally 
that it results not from nuclear effects in the initial state (in particular, gluon saturation) but 
from the final state interaction of hard scattered partons or their fragmentation products in the 
dense medium generated in Au + Au collisions |2()8^ 1209^ \'21{)\ Elj . 

Figure EH shows correlations of high p^ hadrons. The left panel shows the azimuthal distri- 
bution of hadrons with Pt>2 GeV/c relative to a trigger hadron with p7^*"^>4 GeV/c. A hadron 
pair drawn from a single jet will generate an enhanced correlation at A0 ~ 0, as observed for 
p + Pi d + Au and Au + Au, with similar correlation strengths and widths. A hadron pair 
drawn from back-to-back dijets will generate an enhanced correlation at A0 ~ vr, as observed 
with somewhat broader width than the near-side correlation peak for p + p and d + Au coUi- 
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sions. However, the back-to-back dihadron correlation is strikingly absent in central Au + Au 
collisions, and uniquely in central Au + Au collisions. If the correlation is indeed the result of 
jet fragmentation, this suppression is also due to the final state interaction of hard scattered 
partons or their fragmentation products in the dense medium generated in Au-|-Au collisions 
[71] . Finally, the right panel shows the finite azimuthal correlation strength of high pt hadrons 
with the orientation of the reaction plane in non-central Au+Au collisions, in analogy to the 
elliptic flow seen at low pt (Sect. I5.3|l j224j . Since the azimuthal orientation of the initially 
scattered hard parton is uncorrelated with that of the reaction plane and the bulk deformation 
of the fireball, this correlation can only arise from final state interactions. 

Figures EZI and |2H1 present compelling evidence that high p-r hadron production in nuclear 
collisions at RHIC is profoundly altered by interactions with the medium created in the collision. 
We will now discuss in some detail the theory of partonic energy loss, which provides a unified 
description of these phenomena and enables them to be applied as unique, penetrating probes 
of the medium. 



6.1 Partonic energy loss and modified fragmentation functions 

In contrast to the QED energy loss of electrons in matter, the QCD energy loss of partons 
cannot be measured directly because partons are not the final, experimentally observed particles. 
Instead, studies of partonic energy loss must exploit the particle distributions within a jet, in 
particular the modification of the fragmentation functions Da^h{z, n'^) which can be directly 
related to the energy loss of the leading parton. 

The first example we will consider is electron-nucleus deep inelastic scattering (DIS) \2'2'2\ 
l221U225j . We study the semi-inclusive process, e{Li) + A{p) — > e{L2) + h{£h) + ^ , where Li and 
L2 are the four-momenta of the incoming and outgoing leptons, and ih is the observed hadron 
momentum. The differential cross section for the semi-inclusive process can be expressed as 

where p = \p~^, 0, 0_l] is the momentum per nucleon in the nucleus, q = L2—L1 = [—Q^/2q~, q~ , 0± 
is the momentum transfer, s = (p + Li)"^ and oem is the electromagnetic (EM) coupling constant. 
Lfj^u is the leptonic tensor, while W^jy is the semi-inclusive hadronic tensor. 

In the collinear factorization approximation to the parton model, the leading-twist contribu- 
tion to the semi-inclusive cross section can be factorized into a product of parton distributions, 
parton fragmentation functions and the partonic cross section. Including all leading log radiative 
corrections, the lowest order contribution {0{a^)) from a single hard j* + q scattering can be 
written as 

dW^ r 

= T.^lJ dxf^{x,f,j)Hj^J{x,p,q)D,^,{z,,f,'); (52) 

1 271 

Hl°J{x,p,q) = -Tr(7-p7^7- (g + xp)7^)^^5(x-Xij), (53) 

where the momentum fraction carried by the hadron is defined as Zh = d-h/q^ and xb = 
Q"^ /2p^q~ is the Bjorken variable, /ij and /i^ are the factorization scales for the initial quark 
distributions /^(x, /if) in a nucleus and the fragmentation functions Dq^h{zh-, jji^)-, respectively. 

The propagating quark in DIS off a nucleus will experience additional scatterings with other 
partons from the nucleus. The rescatterings induce additional gluon radiation and cause the 
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leading quark to lose energy. This effectively gives rise to additional terms in the evolution 
equation, leading to the modification of the fragmentation functions in a medium. These are 
called higher-twist corrections since they involve higher-twist parton matrix elements and are 
power-suppressed. We will consider those contributions that involve two-parton correlations 
from two different nucleons inside the nucleus. Generalized factorization is usually applied to 
these multiple scattering processes |226t I^TFj I228j . In this approximation, the radiative correction 
to the semi-inclusive tensor from double quark-gluon scattering is 

where 

are twist-four parton matrix elements of the nucleus. The fractional momentum Xi^ = i'^/2p^q^ z{l- 
z) and X = xb = Q'^/2p^q^ is the Bjorken scaling variable. The dipole-like structure in the 
effective twist-four parton matrix results from Landau- Pomeranchuk-Migdal (LPM) interference 
in gluon bremsstrahlung |214l I2l5] . After expansion, the first diagonal term corresponds to the 
so-called hard-soft process where gluon radiation is induced by the hard scattering between the 
virtual photon and a quark at momentum fraction x. The quark is knocked off-shell by the 
virtual photon, returning on-shell by radiating a gluon. The on-shell quark or radiated gluon 
will then have a secondary scattering with another soft gluon from the nucleus. The second 
diagonal term is due to the double hard process where the quark is on-shell after the first hard 
scattering with the virtual photon. The gluon radiation is then induced by the scattering of 
the quark with another gluon that carries finite momentum fraction x^. The two off-diagonal 
terms represent interference between the hard-soft and double hard processes. In the limit of 
coUinear radiation {xl — > 0) or when the formation time of the gluon radiation, tj = 1/xlp~^, is 
much larger than the nuclear size, the two processes interfere destructively, leading to the LPM 
interference effect. 

Including the virtual corrections and the single scattering contribution, we rewrite the semi- 
inclusive tensor in a factorized form with a nuclear modified fragmentation function. 

Ml a, dz 



z 



X 



^lq^qgiz,X,XLjT)Dg^hiZh/z) + /^-fq^gq^Z , X , X lAt)D g^h{,Zh / z) ,(56) 



where Dg^h{zh, f^"^) and Dg^h{zh, fi"^) are the leading-twist fragmentation functions. The modi- 
fied splitting functions are 



'^Iq^ggiz, X, X]^, irp^ 



1 + Z^ 
(1^ 



-TfJx,XL) + S{l-z)AT^(x,fr) 



27iasCA 



A7, 



X, X ij 



5 '"Tl 



A7, 



1 - Z,X,Xl 



, 'i.rpj 



(57) 
(58) 



This medium correction is very similar in form to that caused by gluon bremsstrahlung in vacuum 
that leads to the DGLAP evolution in Eq. Q. 
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Using the factorization approximation |227l I229j . we can relate the twist-four parton matrix 
elements of the nucleus to the twist-two parton distributions of nucleons and the nucleus, 

T,'l{x,xr.) = — (l-e--i/-?^)ff^(a: + Xi)xT/f(xT) + /,^(x)(x^ + XT)/f(xz. + XT) 



Xa 
C 



— {l-e-</'^-)f^{x), (59) 
Xa 

where C is a constant, Xt =< > /2p~^q~z is related to the intrinsic transverse momentum of 
gluons inside the nucleus, xa = ^/rriNRA, fq{x) is the quark distribution inside a nucleus, and 
fg{x) is the gluon distribution inside a nucleon. The coefficient C = 2CxTf^ (xt) should in 
principle depend on and xt but can be approximated as a constant. A Gaussian distribution 
in hght-cone coordinates is assumed for the nuclear distribution, p{y~) = noexp{y~ /2i?^ ), 
where = \/2RAmj\f/p^ and mj^j is the nucleon mass. We should emphasize that the twist- 
four matrix element is proportional to 1/xa = Ra^^n-, ^-e. the nuclear size[229j. 

In the above matrix element, l/xip'^ = 2q~z{l — z)/ ix^ is identified as the formation time 
of the emitted gluons. For formation time that is large relative to the nuclear size the above 
matrix element vanishes, exhibiting the typical LPM interference effect. This results because the 
emitted gluon (with long formation time) and the leading quark remain a coherent system while 
propagating through the nucleus. Additional scattering will not induce more gluon radiation. 

The reduction due to LPM interference of the phase space available for gluon radiation is 
critical for applying the LQS formalism (Luo, Qiu and Sterman j22Hl 1223 122HI) to the problem 
under consideration. In the original LQS approach, the generalized factorization for processes 
with large final transverse momentum £^ ~ leads to consideration of the leading contribution 
in which is enhanced by the nuclear size Ra ~ A^/^. For large Q"^ and A, the higher- 

twist contribution from double parton rescattering that is proportional to agRA/Q"^ will then 
be the leading nuclear correction. Contributions from more than two parton rescattering can 
be neglected. In deriving the modified fragmentation functions, we however have to take the 
leading logarithmic approximation in the limit £y ^ Q^, where It is the transverse momentum 
of the radiated gluon. Since the LPM interference suppresses gluon radiation whose formation 
time (tj ~ Q^/^tP"*") is larger than the nuclear size tjinRa/P^ in our chosen frame, should 
then have a minimum value of i'^ ~ Q'^ /ttinRa ~ Q'^/A^^^, where ttin is the nucleon mass. 
Therefore, the leading higher-twist contribution proportional to Os-Ra/^t ~ otsR\/Q'^ from 
double scattering depends quadratically on the nuclear size Ra- 

With the assumption of the factorized form of the twist-4 nuclear parton matrices, there 
is only one free parameter C(Q^) which represents the quark-gluon correlation strength inside 
nuclei. Once it is fixed, the energy and nuclear dependence of the medium modification 
of the fragmentation function can be predicted. Shown in Fig. 1221 are the calculated nuclear 
modification factors of the fragmentation functions for ^^A^ and ^'^Kr targets, compared to recent 
HERMES data P^ l23Uj . The predicted shapes of the z- and //-dependence agree well with the 
data. A remarkable feature of the prediction is the quadratic A?/'^ nuclear size dependence, 
which is verified for the first time by an experiment. By fitting the overall suppression for one 
nuclear target, one fixes the only free parameter in the calculation, C{Q'^) = 0.0060 GeV^ with 
as(Q^) = 0.33 at g2 ~ 3 GeV^. 

Modification of the fragmentation can be quantified by the quark energy loss, defined as the 
momentum fraction carried by the radiated gluon: 

{AZg){xB,^l'^) = / ^ / dz'^zA-fq^gq{z,XB,XL,iT) 

Jo It Jo Zn 
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Figure 29: Measured ratios of hadron distributions from DIS off A and d targets, from HERMES 
|23mi23l] . The curves show the predicted nuclear modification of the jet fragmentation function, 
described in text. Left: vs. fragmentation fraction z. Right: vs. energy transfer v. 



Nr 



Jo It jo 



T^g{xB,XL) 



C 



CaoI xb 

iVc XaQ"' 



'I 



dz^——^ 

z{l — z 



z 



dxi 
xl 



(60) 
(61) 



where x^ = jjp' /2p^q z{l — z) = xb/z{1 — z) for factorization scale /i^ = Q^. For xa ^ xb ^ 1, 
the leading quark energy loss is roughly 



(A2;^)(xB,/i') ~ C 



CaoI Xb 



6v^ln 



2XB 



(62) 



Since xa = ^/m^RA, the energy loss (A^;^) thus depends quadratically on the nuclear size. 

In the rest frame of the nucleus, = ttin, q~ = z/, and xb = Q^/2p+g~ = Q^/2mNt^- 
The average total energy loss is AE = u{Azg) ^ C{Q^)a'^{Q'^)mNRAlCA/Nc)3y/7ilia{l/2xB). 
With the value of C from the fit, (xb) ^ 0.124 in the HERMES kinematics [2201 ESI] and the 



average distance {La) = RAy^/ir for the assumed Gaussian nuclear distribution, the average 
quark energy loss dE/dL ^ 0.5 GeV/fm in a Au nucleus. 

Attenuation of leading hadrons in DIS off nuclear targets has also been studied in hadronic 
transport and absorption models |232| I233j . In these models, two distinct types of hadronic 
absorption are assumed: absorption of fully formed physical hadrons, and of "pre-hadrons" . 
Since the hadronic formation time is long, attenuation effects on fully formed hadrons are small 
and the observed attenuation is attributed to the absorption of pre-hadrons in the nuclear 
medium. The pre-hadron can be modeled as a quark-antiquark dipole, and the interaction of qq 
dipoles with the nuclear medium should be equivalent to the picture of multiple parton scattering 
and induced bremsstrahlung |220j . 



6.2 Energy loss and jet quenching in a hot medium 

To extend the study of modified fragmentation functions to jets in heavy-ion collisions, we 
assume a one-dimensional boost invariant (Bjorken) expansion with transverse gluon density 
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profile p{r,r) = {tq/t)6{Ra — r)po, and (k^) ~ (the Debye screening mass). The initial jet 
production rate is independent of the final gluon density, which is related to the parton-gluon 
scattering cross section asXTG^Xx) ~ p^Cg |234j so that 



XL 



xb 



I? j dragp{r,T)[l-cos{r/Tf)], 



(63) 



where Tf = 2Ez{l — z)/l\ is the gluon formation time. Assuming partons traveling in the 
transverse direction at the velocity of light (r = r — tq), the fractional energy loss from Eq. ()61|) 
is L235i 
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Keeping only the dominant contribution and assuming ag 
gg), the average energy loss is 



Ca2TTal/p^ {Ca=l for qg, 9/4 for 



dL Ra Jto 



dTp{T){T - To) In 



2E 



(65) 



Neglecting the logarithmic dependence on r, the averaged energy loss in a one- dimensional 
expanding system is {dE/dL)id ~ {dEQ/dL){2To/ Ra), where dE^/dL oc PqRa is the energy loss 
in a static medium with the initial gluon density po of the expanding system at time Tq. Because 
of the expansion, the averaged energy loss {dE/dL)id is suppressed relative to the static case 
and does not depend linearly on system size for a fixed value of po- 

This form of the energy loss has also been derived in the opacity expansion framework 
|236| 123 7j. based on a model of multiple scattering in a quark-gluon plasma consisting of ran- 
domly distributed scattering centers with screened static potential [212J. The magnitude of the 
momentum transfer is small, limited by the Debye screening mass p ~ gT, and amplitudes for 
multiple scattering and gluon bremsstrahlung factorize in momentum space. This leads to an al- 
gebraic reaction operator formulation of the radiation amplitude induced by multiple scattering, 
which keeps track of the phase accumulation due to multiple scattering. This approach enables 
the iterative evaluation of the radiation spectrum induced by a given number of scatterings, 
corresponding to an expansion in opacity parameter ^ = apL. However, this approach cannot 
be applied directly to multiple parton scattering in a cold nuclear medium since the concept of 
random screened potential is difficult to justify in that case, and multi-parton correlations are 
also important. The extension of the twist expansion to higher twist to account for multiple 
parton scattering is also difficult and has yet to be done. To first order in opacity or leading 
twist, the two approaches give the same result for radiative energy loss in a hot gluon plasma. 

In order to calculate the effects of partonic energy loss on the production of high px hadrons 
in nuclear collisions, one can apply a simpler effective modified fragmentation function |238t I231J] . 



D'^^^iz,,Q^AE,) = il-e 



z, 



A Zc 

+ e~^^^Dl/^{z.,Q% 



(66) 



where z'^, Zg are the rescaled momentum fractions. The fragmentation functions in free space 
^h/ci^cyQ"^) are given by the BBK parameterization |24Uj . The first term is the fragmentation 



52 



function of the jet c after losing energy AEc{pc, 0) due to medium induced gluon radiation. The 
second term is the feedback due to the fragmentation of the Ng{pc, (p) = {AL/X) radiated gluons. 
This effective model is found to reproduce the pQCD result from Eg . (j56j) very well, but only 
when Az = AE^/E is set to be Az ~ 0.6(2g). Therefore the actual averaged parton energy loss 
should be AEjE = I.6A2; with Az extracted from the effective model. The factor 1.6 is mainly 
due to the effect of unitarity correction in the pQCD calculation. 

Since gluons are bosons, there should also be stimulated gluon emission and absorption by the 
propagating parton because of the presence of thermal gluons in the hot medium. The detailed 
balance is crucial for parton thermalization and should also be important for calculating the 
energy loss of an energetic parton in a hot medium. Taking into account such detailed balance 
in gluon emission, the asymptotic behavior of the effective energy loss in the opacity expansion 
framework is |241j . 



absorption in detailed balance, which effectively reduces the total partonic energy loss in the 
medium. 

Fig. inOl shows numerical calculations of the ratio of radiative energy loss with and without 
stimulated emission and thermal absorption, as a function of E/^ for L/Xg = 3,5 and = 0.3. 
The insert shows the energy gain via gluon absorption with and without rescattering. For very 
high energy partons the effect of the gluon absorption is negligible. However, for intermediate 
parton energy the thermal absorption reduces the effective parton energy loss by about 30-10%, 
generating an energy dependence of the effective parton energy loss in the intermediate energy 
region. This energy dependence is parameterized as 



The threshold is the consequence of gluon absorption, which competes with radiation to effec- 
tively shut off the energy loss. According to Eq. ()65|). the parton energy loss and therefore the 
parameter eo in the above equation depends linearly on the gluon density of the medium. This is 
the only property of the medium one can extract from the experimental measurement of parton 
energy loss. 

6.3 Cronin Enhancement in j9 + ^4 Collisions 

Tomography relies upon accurate knowledge of the initial flux of radiation in order to measure the 
density of the intervening matter between the source and the detector. Jet quenching is suitable 
for tomography since jet production can be calculated reliably in pQCD, in good agreement 
with experimental measurements in high-energy p + p{p) collisions. Perturbative calculations 
and experimental measurements of jet production cross sections in p + p and p + A collisions 
effectively calibrate the initial source of the beam of jets. 

The simplest observable of jet production in p + p collisions is the high pt hadron spectrum 
resulting from jet fragmentation. In the coUinear factorized parton model, the single inclusive 
hadron spectrum is |242j 





i-rrhd = eo(i^//i - 1.6)1-7(7.5 + E//i), 






■PT 



I dxadxbd'^kaTd^kbTgp{kaT,Q'^)gp{hT,Q'^) 




53 




Figure 30: The ratio of effective parton energy loss with absorption {AE = AE^^^ + AE^j^J + 
AE^^J^) to that without (AE^^J^), as a function of E/n. Insert: energy gain via absorption with 
rescattering (Ai?^),^) and without (Ai?^°^). 



X fa/p{Xa,Q'^)fb/p{Xb,Q 



TlZr 



dt 



{ah — > C(i), 



(69) 



where fa/p{x) is the parton distribution in the proton (we will use MRSD-/ parameterization 
|248j ) ■ D^^^{zc,Q'^) is the fragmentation function of parton c into hadron h derived from e~^e~ 
data |24()j . and Zc = Pt/ptc is the momentum fraction of the jet carried by a produced hadron. 
The K 1 — 2 factor accounts for higher order QCD corrections. At fixed-target energies, NLO 
calculations with resummation underestimate the measured hadron spectra |244j . indicating the 
importance of higher power corrections. This can be remedied phenomenologically by introduc- 
ing an intrinsic transverse momentum smearing 5'p(A;r,Q^), assumed to have a Gaussian form 
dpikx) = e~^r/(^T)/7r(/c|,). This smearing is found to be important for < 100 GeV [89j. It 
is less important at collider energies as seen in Fig. |^ which shows good agreement between 
LO and NLO parton model calculations without kx and measured charged hadron and vr'' single 
particle inclusive spectra. Additional discussion of hadronic spectra and intrinsic kx can be 
found in piKlIM)] . 

Two nuclear effects must be incorporated into the parton model for an accurate description 
of hadron production in p + A collisions, both arising as a consequence of multiple scattering: 
the nuclear modification of the parton distribution functions and nuclear kx broadening. These 
effects must also be taken into account for quantitative study of the change in hadron production 
due to jet quenching in A + A collisions. 

A practical approach to multiple scattering effects calculates the differential cross section for 
inclusive hadron production as that for a single hard parton-parton scattering, but with larger 
beam parton kx due to multiple soft partonic scattering prior to the hard interaction and parton 
distributions modified according to DIS measurements off nuclear targets. The single inclusive 
hadron production cross section in minimum-bias p + A collisions is then 
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Figure 31: Single particle inclusive spectra at collider energies compared to parton model and 
NLO calculations. Left: charged hadrons from p + p collisions at y/s = 200, 900, 1800 GeV from 
UAl |247j and CDF |248j . Lines indicate parton model calculations with (solid) and without 
(dot-dashed) intrinsic kr- Right: vr^ from p + p collisions at y^=200 GeV compared to NLO 
calculations, from PHENIX |249j. 

TTZc at 

where TA{b) is the nuclear thickness function normalized to / (PbTA{b) = A. The parton distri- 
bution per nucleon inside the nucleus of mass A and charge Z at impact parameter b, 



' Z Z 

fa/A{x, Q'\ b) = Sa/A{x, b) -jfa/pix, Q^) + (1 - ■j)fa/n{x, Q^) 



(71) 



is assumed to be factorizable into the parton distribution in a nucleon fa/N{x, Q"^) and the nuclear 
modification factor Sa/A{x, b), parameterized in various ways |250t I25l[ll02j according to the DIS 
data. 

The initial partonic transverse momentum distribution in a projectile nucleon striking the 
target nucleus at impact parameter b is still assumed to be Gaussian but with a broadened 
variance 

(4)a(Q') = + S'miMb) - 1). (72) 

The broadening is assumed to be proportional to the mean number of scatterings z/a(&) the 
projectile suffers inside the nucleus. The parameters are fitted to existing fixed-target p + A 
data at energies up to i/i = 40 GeV [5^ 

For quantitative comparison of hadronic spectra from A + B and p + p collisions we define 
the nuclear modification factor 



Rab{pt) = . , „ ■ (73) 

{Nun)daNN dyd^PT 
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Figure 32: Rab{pt) (eq. ITHj) for midrapidity vr" and charged hadrons from d+Au relative to 
p+p collisions at 200 GeV, from PHENIX [203 and STAR [H]. Left: data compared to first 
predictions of the Cronin effect in p + Au collisions at y/s = 200 GeV [HH]- Right: data compared 
to calculated Cronin effect for 7r° from a Glauber-eikonal model of multiple scattering |253j . The 
solid and dashed lines show variation with intrinsic kr, dotted band shows uncertainty due to 
infrared regulator pq = 1.0 ± 0.1. 



{Nhin), the mean number of binary collisions, is discussed in Section 13.31 In the absence of 
nuclear effects the cross section for a hard process will scale as N^in, making Rab{pt) unity. 
Nuclear-specific effects are measured by the deviation of Rab{pt) from unity. In the following 
we will also utilize a related quantity Rcp{pt), the binary-scaled ratio of central over peripheral 
spectra from A-fA collisions: 

R^pipj,) - ^^^/ dPTdrj/ (Nhin) (central) 
cPN/ dpTcLr]/ {Nbin) (peripheral) 

Rab(pt) and Rcp(pt) contain similar physics, though Rcp(pt) niay at times be preferable from 
the standpoint of experimental uncertainties. 

Fig. 1221 shows Rab(pt) for 7r° and charged hadrons from d-|-Au collisions at 200 GeV, mea- 
sured by PHENIX [209j and STAR |74j (see also PHOBOS EJU] and BRAHMS [201 )• The 
curves in the left panel show the first prediction of the nuclear modification of hadron spectra in 
p + A collisions within the parton model p3j. The enhancement at intermediate px, known as 
the Cronin effect, is due in this calculation to broadening of kj- from multiple scattering. The 
enhancement disappears at large px, as do all higher-twist effects. The modest Cronin enhance- 
ment predicted for p + A collisions at this energy is confirmed by the measurements. The Cronin 
effect in A -|- A collisions at this energy should be of similar magnitude. 

The Cronin effect has also been studied in multiple parton scattering models which utilize 
the eikonal Glauber framework [206| I252| I253j . Interference effects play an important role, 
for instance the absorption arising from single scattering cancels part of the double scattering 
contribution, leading to a 1/pt^ dependence of the yield enhancement due to double scattering 
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Figure 33: High pt single inclusive hadron spectra from centrality-selected Au+Au collisions at 
^s;;;;=200 GeV. Left: charged hadrons from STAR [S^. Right: 7r° from PHENIX 



|2()6j . The dominant double scattering contribution occurs when one scattering is soft and most 
of the jet's pt comes from one hard scattering. At lower pt the absorptive correction dominates 
double scattering and the coherence of the interaction suppresses the spectrum relative to binary 
scaling. The transition between the coherent suppression and Cronin enhancement occurs at a 
scale po, where parton-nucleon scattering is no longer a power-law like process. Fig. 1221 right 
panel, shows such a calculation |253j compared to data. 



6.4 High Pt hadron suppression in ^4 + ^4 coUisions 

The foregoing analysis of hard scattering and high pt hadron production at midrapidity in p + p 
and p{d) + A collisions demonstrates that the elementary production cross sections and nuclear 
multiple scattering effects are well understood, providing a foundation for jet tomographic studies 
in A + y4 collisions. Within the energy loss picture of multiple parton scattering and induced 
gluon bremsstrahlung, high pt jet fragmentation will be softened and the final hadron spectra 
suppressed in A + A collisions |2U2j . The effective partonic energy loss can therefore be deduced 
from measurements of high pt hadron suppression, providing a direct measurement of the initial 
gluon density. The experimental input for this study is shown in Fig. EHl Charged hadron and 7r° 
single particle inclusive spectra from centrality-selected Au + Au collisions have been measured 
with high precision over a very broad Pt range |254L I255j . 

The inclusive hadron spectrum in ^4 + A collisions can be calculated using a LO pQCD model 



da f 
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X fa/Aixa, Q ,r)fb/Aixb,Q ,\h-r\) — -^{ah^cd), (75) 

nzc dt 

with medium modified fragmentation functions -D^/^ given by Eq. 
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Figure 34: Binary collision-scaled ratios of spectra shown in Fig. in^|254[ I255j fsee also PHOBOS 
j257l l2Tn] and BRAHMS ^Q^). Left: Centrality dependence of Raa{pt) vs. pt- Theory 
curves are from calculations described in text. Boxes at Raa{pt)=^ indicate the (multiplicative) 
uncertainty due to binary collision scaling. The hatched region on the central collision theory 
curve indicates variation in eo of ±0.3 GeV/fm (Eq. lUHj) . Right: Rcp{pt) [Eq- (IZH)] for central 
(0-5%) over peripheral (40-60%, 60-80%) spectra, for charged hadrons from STAR |254j . 



Assuming one-dimensional longitudinal expansion and gluon density Pg{T, r) proportional to 
the transverse density of participant nucleons, the impact parameter and path length dependence 
of the energy loss is [Eq. (jHSI)] 

AE(6,r,0) ^ (— )i, / dr '-p^(r,b,r + nr), (76) 

dL Jto tqPo 

where AL(6, r, 0) is the path length in matter for a jet produced at r and traveling in direction 
n with azimuthal angle relative to the reaction plane, in a collision with impact-parameter b. 
Po = {Pg)ijo) is the average initial gluon density in central collisions at tq and {dE/dL)id is the 
average parton energy loss in a one-dimensionally expanding medium with an initial uniform 
gluon density po- The corresponding energy loss in a static medium with a uniform gluon density 
Po within radius Ra is |235j dE^/dL = {RA/2TQ){dE/dL)iii. We will use the parameterization 
in Eq. for the effective energy dependence. 

Fig. OH left panel, shows Raa{pt) derived from the data in Fig. 1221 with additional 
charged hadron measurements from PHENIX ^258j . The figure shows comparable charged hadron 
Raa{pt) from STAR (filled circles) and PHENIX (triangles), which are seen to differ in overall 
normalization by ~20%. This difference can be attributed to the ~20% difference in normaliza- 
tion of the p + p reference spectrum used in the two analyses (denominator of Eq. ITHjl |254| 1258] . 
The measurements are however consistent within systematic uncertainties and for the purpose of 
this discussion the difference can be taken as an indication of the precision of the measurements. 
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The figure shows a marked effect: high pt hadron production in central Au + Au collisions is sup- 
pressed by a factor five relative to binary collision scaling. The suppression has strong centrality 
dependence, with hadron yield in the most peripheral collisions consistent with binary collision 
scaling. Similar results have been obtained by PHOBOS [2^31^01 and BRAHMS [203 • This 
suppression stands in strong contrast to the enhancement in high px hadron production observed 
for d + Au collisions in Fig. |2ZI and 1221 The strong yield suppression in central Au+Au con- 
trasted with the enhancement in d+Au collisions clearly demonstrates that the central Au+Au 
suppression at high pt is due to final-state interactions in the fireball generated by the collision. 

The solid lines in Fig. |^ left panel, are parton model calculations according to Eq. ()75p. 
The fit to the observed factor of five suppression at large pt in the more central collisions results 
in the parameters eo = 1.07 GeV/fm, /i = 1.5 GeV for the effective quark energy loss in Eq. (|UHj). 
The difference between charged hadron and 7r° spectra at intermediate pt is attributed to a non- 
perturbative component of hadron production, which we will return to in the next subsection. 
To demonstrate the sensitivity to the parameterized partonic energy loss in the intermediate 
Pt region, the dashed curves show Raa{pt) for charged hadrons excluding the soft component 
|256j . Once these parameters are fixed by fitting to central collisions, the centrality dependence 
of the suppression is a prediction of the model. 

The suppression factor at the highest measured pt is shown for charged hadrons in Fig. 
IMl right panel |254j {Rcp{pt) extends the pt reach and reduces the systematic uncertainties 
relative to Raa{pt) |254j ). The suppression again is seen to be a factor five, independent of 
Pt- This in) dependence is well reproduced by the pQCD calculations incorporating partonic 
energy loss (pQCD-I [256. , pQCD-II [259j). The pT-independence of the suppression is not a 
trivial effect: it results from the subtle interplay between Cronin enhancement, shadowing of the 
structure functions, and jet quenching in a parton model with energy loss |259j . 

The suppression for vr" (left panel) varies only slightly over a very broad pt range, down 
to Pt ~1-2 GeV/c. A constant or logarithmic jet energy dependence of the energy loss as in 
cold nuclear matter generates a suppression factor that slowly rises with pt |l235j . Thus, the 
energy dependence of the partonic energy loss seen in A+A collisions is quite different from that 
in DIS off cold nuclei. This is an indication of detailed balance at play in a thermal system. 
Detailed balance also leads to a slight rise in the calculation of Raa{pt) at pr<4 GeV/c, where 
the fragmentation picture gradually loses its validity and hadron production may be influenced 
by non-perturbative effects, especially for kaons and baryons. For pt>5 GeV/c, these non- 
perturbative effects are expected to diminish and indeed the tt" and charged hadron suppression 
factors converge in that region. 

Partonic energy loss calculations reproduce the centrality dependence of the suppression 
well. The solid curves in the left panel of Fig. |^ are calculated using partonic energy loss 
from Eq. (jUH)) with parameters fit to the central collision data for vr" and charged hadrons. The 
centrality dependence is determined by the averaged total energy loss [Eq. (f7B)) ]. which results 
in an effective surface emission of the surviving jets. Jets produced at the core of the overlap 
region are strongly suppressed, since they lose the largest amount of energy. The centrality 
dependence of the suppression is found to be dominated by the geometry of the produced dense 
matter rather than the length dependence of the parton energy loss. The non-Abelian nature 
of the path length dependence of the energy loss must be tested in other ways, in particular via 
central A + A collisions with varying nuclear size. 

Transverse expansion of the bulk medium can also be considered in the parton model calcu- 
lation. Since the total energy loss in Eq. ()65|) depends on a path integral of the parton density 
profile p(r, r), the transverse expansion will increase the duration of parton propagation in mat- 
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Figure 35: Binary collision-scaled ratio of yields Rcp{pt) (eq. 17^ from central relative to 
peripheral collisions for identified particles. Left: protons and pions from PHENIX |76j. Right: 
strange mesons and baryons from STAR I264j compared to charged hadrons from Fig. |23 



ter. However, the expansion will also accelerate the reduction of the parton density along the 
path. These effects compete and the final total parton energy loss remains approximately the 
same as in the case without transverse expansion |2tiUj . 

It was proposed recently that gluon saturation effects (Sect. 14. ip can extend well beyond 
the saturation momentum scale Qs, resulting in hadron suppression relative to binary scaling 
{Rab{pt)<^) for hadron p-p ~ 5 — 10 GeV/c at RHIC energies j261j as shown by the solid lines 
in the right panel of Fig. |^ Since this suppression originates in the properties of the incoming 
nuclear wave-function, hadron production in d + Au collisions should also be suppressed by this 
mechanism |2tilj . Experimentally, however, an enhancement in mid-rapidity hadron production 
md+Au is seen instead (Figs. l^and lHipiinH7Ill^lM?j ). even in central d+Au collisions jH] 
where saturation effects should be most pronounced. The observed enhancement is at variance 
with saturation model expectations |261j . 

The saturation model calculation has since been developed to include the Cronin effect due to 
classical elastic scattering |262j . Quantum evolution cancels the Cronin enhancement, however, 
leading to suppression of high pt hadrons at high energies or large rapidities. Suppression at 
large rapidity might however not be a unique signature of gluon saturation. For instance, limiting 
fragmentation of valence quark jets could also result in hadron suppression at forward rapidity 
|2ti3j . Discrimination of these mechanisms may be possible by utilizing di-hadron correlations, 
which we discuss below. 



6.5 Parton Recombination 

The hadron suppression factor Raa{pt) in Fig. |^ is different for charged hadrons and in 
the intermediate pr ~ 2 — 5 GeV region. More detail of this effect is seen in Fig. EH The 
left panel shows that 7r° are strongly suppressed in central collisions whereas protons are not 
suppressed in this region. The right panel shows similar systematic behavior for strange meson 
and baryon yields. It will also be seen that the azimuthal anisotropy (^2) depends on hadron 
flavor at intermediate pt-, with the anisotropy of baryons larger than that of mesons |14tij . 
In the intermediate pt region, where the parton energy loss is strongly influenced by the 
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medium through absorption of thermal partons, fragmentation or hadronization should also be 
modified by the presence of other partons. These nonperturbative effects invalidate the picture of 
independent parton fragmentation in the vacuum. For instance, the leading parton may pick up 
another parton from the medium before hadronization, a scenario referred to as recombination 
or coalescence. Such recombination processes were long ago proposed as the dominant mech- 
anism for hadron production in the beam fragmentation region of hadron collisions ^265| I266j . 
Since recombination is a non-perturbative process, hadron wavefunctions should be taken into 
account along with the initial parton distributions. Alternatively, the recombination processes 
can be described as higher-twist corrections to the fragmentation functions, since they involve 
higher-twist matrix elements of the overlap between two- or three-parton operators and the fi- 
nal hadronic state. Such an approach has been applied to D meson production to explain the 
— D~ asymmetry in the forward rapidity region oi h + p collisions |267 |. 




To take this effect into account in the simple parton model a nuclear dependent soft compo- 
nent, assumed to be proportional to (A^binary), is added to kaon and baryon fragmentation func- 
tions so that in central Au + Au collisions the ratio {K +p)/n 2 at ~ 3 GeV/c and the ratio 
approaches its value ioi p + p collisions at pt > 5 GeV/c. The resulting suppression for charged 
hadrons and its centrality dependence agree well with the STAR data (Fig. IMj) . The and 7r° 
suppression are related via the (A' + p)/7r ratio: i?^^ = R\j^[l + {K + p)/n]AA/['^ + {K + p)/7!-]pp. 
It is apparent from the data that {K + p)/ti converges for Au + Au and p + p colhsions at px > 5 
GeV/c. Since such nonperturbative effects are caused by presence of other produced partons, 
thermalized or not, a qualitatively similar dependence of the Cronin enhancement on hadron 
species should be seen in p{d) + A collisions. Such an effect is indeed observed |268J, though the 
enhancement relative to p + p collisions in the ratio of proton to pion yields at intermediate pt 
is markedly smaller in ci -|- Au than in Au + Au collisions. 

The observation of the flavor dependence of the suppression factor and, in particular, the 
flavor dependence of the azimuthal anisotropy (Sect. 16. 7|) have spurred new developments in 
the recombination approach to hadron production in heavy-ion collisions (269^ . Parton recom- 
bination in heavy-ion collisions was first investigated by Hwa and Yang |27m I271j and was later 
studied in detail in terms of coalescence or recombination models |272t 12731 1274| I275j . These 
models generally assume two components of high pt hadron production in heavy-ion collisions. 
The recombination of partons from the bulk medium dominates the production of low and in- 
termediate Pt hadrons, while high pt hadrons result mainly from the fragmentation of parton 
jets after propagating through the bulk medium and losing energy. In these models, the number 
of mesons formed from parton recombination is |275j 



where da denotes the differential element of the space-like hypersurface of hadronization and 
fM{xi,X2;pi,P2) is the coalescence probability given by the Wigner distribution function of the 
meson in terms of constituent quarks. The statistical factor qm takes into account the number 
of internal quantum states in forming a colorless meson from a colored quark and antiquark. 
The formula for baryon production from parton recombination is similar, except that it involves 
the baryon coalescence probability from three constituent quarks. The quark (or antiquark) 
distributions should contain both soft and hard components. The soft component arises from 
a thermalized quark-gluon plasma with collective radial and elliptic flow though the connection 
between the consituent quarks appearing in Eq. (f77j) and the massles partons of a chirally 
restored plasma is not specifled at present in the model. The hard component is given by the 




Nm = 9m / Pi ■ daip2 ■ da2 



d^Pi d^p2 



fq{Xl]Pl)fq{X2] P2)fM{Xi, ^1, ^2), 



(77) 
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Figure 36: pt dependence of the ratio of baryon to meson yields at ys^=200 GeV. Left: p/n~ 
for central Au + Au collisions, from PHENIX [TB] . Lines are recombination model calculations 
[276^ with (solid) and without (dashed) soft+hard contributions due to coalescence of minijet 
partons with thermalized partons. Right: A/K^ for p + p and centrality-selected Au + Au 
collisions, from STAR 



minijet distribution from a pQCD calculation, with transverse momentum reduced in accordance 
with the observed partonic energy loss. 

In the low and intermediate pt regions, the hadron spectra are dominated by recombination 
of thermal partons. Such a description may be considered as a model for hadronization of 
the quark gluon plasma. In this region the hadron spectra are determined by the underlying 
parton spectra and the statistical factors Qm^b- If an exponential form is assumed for the 
parton spectra, the proton to pion ratio in a simple recombination model will be independent 
of transverse momentum and is determined only by the ratio of statistical factors QbI Qm |274j . 
If the contribution of A-decay to the proton yield is included and the contribution of higher 
resonance decay to pions is excluded, p/vr~ ~ 1. Fig. EHl left panel, shows the measured p/vr~ 
ratio for central Au + Au collisions at y^s^=200 GeV from PHENIX [70] compared to the 
calculation by Greco, Ko and Levai |276j . one of several recombination models |270L I274j . The 
calculation describes the data well. The p/vr ratio indeed reaches a value of about unity, much 
larger than the value achieved m p + p collisions j2(j8J and in gluon jets from e^e~ annihilation 
j277j . The pt dependence at lower pt can be attributed to the contribution to the pion spectra 
from resonance decays, whose importance diminishes at higher pt- Recombination models also 
successfully describe azimuthal asymmetry measurements, as discussed in Sect. 16.71 

At large transverse momentum the underlying parton spectrum is dominated by pQCD 
minijet production, which has power-law pt dependence. The recombination mechanism in this 
case will always generate higher meson than baryon yields. For increasing pt-, fragmentation will 
eventually dominate recombination and the baryon/meson ratio should revert to that of parton 
fragmentation in vacuum, consistent with the expectation that recombination is a higher-twist 
process that is suppressed at high p^. This p^-dependence of the relative yield of mesons 
and baryons is common to a variety of recombination models |27m I274[ I275j . It is in broad 
agreement with the data in Fig. EHl right panel, which shows a large enhancement of A/K^ for 
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more central Au + Au collisions at intermediate pt but convergence of the A/K^ ratio for all 
Au+Au centralities and p+p collisions at higher pt- 

The model of Greco, Ko and Levai includes an additional pick-up process by leading hard 
partons, corresponding to the recombination of a parton from the hard component and a parton 
from the soft component |275j . Such soft+hard recombination can be significant at intermediate 
Pt as shown by the dashed line in Fig. IHEl left panel. This process will be critical for addressing 
the observed two-particle correlation of the same-side jet that will be discussed in the next 
subsection. 

Recently, a jet fragmentation model was proposed |278j in which parton recombination is the 
hadronization mechanism even in vacuum. The medium effect at intermediate pt is a natural 
consequence of the hadronization in medium in this model, which also leads to flavor dependence 
of the Cronin enhancement in d + Au collisions [2721 l28Uj . Quantitative tests of recombination 
models and the determination of their parameters are best done using dihadron correlations. 
These studies are now underway. 

As discussed in previous sections, measurements of low pt inclusive spectra and elliptic flow 
can be reproduced well by calculations incorporating ideal hydrodynamics. It is to be expected 
that such a hydrodynamic picture will break down at high pt where partonic cross sections are 
small and the mean free path is long, meaning that hard partons cannot be brought to thermal 
equilibrium as readily as soft partons. In this region, pQCD dynamics should dominate. To 
study the transition between these regimes at intermediate px, hybrid calculations have been 
carried out that combine hydrodynamics with a pQCD parton model |281j . Since collective flow 
gives higher mass hadrons a larger px boost, it is expected that the effects of collective flow 
will extend to larger p^ for baryons than for light mesons. Such a two-component approach 
can indeed describe the observed p/n and A/K^ enhancements in Au + Au collisions. However, 
since the effects are purely kinematic it also predicts the same enhancement for heavy mesons. 
Preliminary data ( |264t I^H^ I283j . see also Fig. EHl right panel) indicate that the enhancement is 
rather more dependent on whether the hadron is a meson or baryon than on its mass. However, 
heavier mesons such as the may have smaller hadronic cross sections than baryons and not 
flow effectively with the bulk medium during the hadronic stage. An interesting test will be the 
measurement of azimuthal anisotropy of (j) meson spectra, which may be relatively insensitive to 
the hadronic dynamics at the late stages of the expansion. 

6.6 Dihadron Spectra and Jet Quenching 

Jets are produced in pairs in leading order pQCD, and high pt dihadron correlations in p + p 
collisions should exhibit the back-to-back jet structure of the underlying hard parton-parton 
scattering. Jet quenching due to partonic energy loss in nuclear collisions is expected to modify 
such correlations. 

As we have shown in Fig. |2H1 left panel, the relative azimuthal angle distributions of high p^ 
hadron pairs in p + p and d + Au collisions exhibit the two-peak feature characteristic of back- 
to-back jet pairs but the away-side correlation vanishes in central Au + Au collisions, consistent 
with the predicted phenomenon of mono-jet production due to jet quenching |284|. In these 
measurements the trigger hadron has pt*"^>4 GeV/c, with the relative azimuthal angle plotted 
for all other hadrons having 2<p'p<pT*"^ |223U74j . All particles lie within |?7|<0.7. For Au + Au 
collisions the effects of elliptic flow on the dihadron distribution must be taken into account 
[223j . The distributions are normalized to the number of trigger hadrons, thereby measuring the 
probability to find an associated hadron. The figure shows background-subtracted distributions. 
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where the background is defined as the yield in an azimuthal interval orthogonal to the trigger 
where the correlated yield in p + p collisions is seen to be small. 

The strength and width of the small-angle peak (A0 ~ 0) are similar in all three distributions. 
The relative probability that the small-angle hadron pair has the same vs. opposite charge sign 
is similar to that measured in jets from + e~ annihilation ("charge ordering" |285| I223j ). 
The contribution of resonances to this peak is estimated to be negligible |223j . This series 
of evidence leads to the conclusion that the small-angle correlation in all systems from p + p 
to central Au + Au results dominantly from the fragmentation of jets. The large-angle peak 
(A(/) ~ tt) has similar strength and width in p + p and d + Au collisions, indicating a hadron 
pair produced by the fragmentation of a back-to-back jet pair. The large-angle high pt dihadron 
correlation is however markedly suppressed in central Au + Au collisions, suggesting significant 
suppression of the back-to-back jet yield. 

The strength of the correlation in Au + Au can be quantified by integrating the yield in the 
correlation peak and comparing to the p + p correlation via |223j 

/ft^ d (A0) -B(l + 2vl cos (2A0))" 

Iaa{pt) = ^ ^ ^ , (78) 

/a^ (A0) C^'^' 

where C^"^" and C'P'p are the measured dihadron azimuthal angle distributions in Au + Au and 
p + p prior to background subtraction. The subtracted background i? (1 + 2f| cos (2A(/))) results 
from a fit to the Au + Au distribution in the background interval orthogonal to the trigger 
direction. The elliptic fiow coefficient f 2 is determined by independent measurements |286j . The 
integration interval [A</)i, A02] spans either the small- angle or back-to-back peaks in the p + p 
distribution. Figure IH7| left panel, shows the centrality dependence of Iaa{pt) for both same-side 
(upper) and back-to-back (lower) dihadron pairs. As is also apparent in Fig. |2H1 the same-side 
correlation strength for all Au + Au centralities is similar to p + p. In contrast, the back-to-back 
correlation strength in Au + Au is seen to vary markedly from peripheral collisions, with strong 
correlation, to central collisions where it is negligible [2.23J. It is compelling to attribute the 
back-to-back suppression in central collisions to jet quenching. 

The population of jets contributing to the near-side and back-to-back correlations may differ 
for two reasons. First, the requirement of a pair of hadrons above threshold biases the near-side 
contribution to jets with larger initial energy. Second, if the partonic energy loss in the core of 
the fireball is large, the trigger will bias towards those jets that are produced near the surface of 
the reaction volume, heading outward. The jets that generate the same-side correlation would 
therefore punch through the medium after losing a limited amount of energy and fragment 
essentially in vacuum, leading to similar near-angle dihadron correlations from p + p to central 
Au + Au, as observed. The jets recoiling against the trigger will however be biased towards the 
population heading into the core of the fireball and will be strongly suppressed, as observed. The 
total energy of the jets, suppressed or not, is of course conserved but it may be redistributed 
to softer hadrons with much broader angular distributions. In the picture of partonic energy 
loss due to induced radiation, the emitted gluons could further interact with other partons and 
dissipate their energy to the medium. The hadronic structure of the jet will thereby be lost, 
leading to strong suppression of the back-to-back high pt dihadron correlation. This picture of 
trigger bias is also consistent with the suppression of inclusive production at high pt relative to 
binary collision scaling (Fig. IH^ . since jets will be emitted dominantly from the surface and not 
the volume of the fireball. 

It should be noted, however, that the correlations discussed thus far have pt*'^'^>4 GeV/c. 
While the systematic features of near-angle correlations in central Au + Au are similar to those 
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Figure 37: High dihadron distributions from 200 GeV Au + Au and p + p collisions, from 
STAR. Left: centrality dependence of near-angle and back-to-back dihadron correlation strength 
in Au + Au relative to p -|- p {Iaa{pt), Eq- [ZH]) for pT{trig)>A GeV/c |223j . Central collisions 
correspond to large Npart- Right: dihadron azimuthal distributions for pT{trig)>Q GeV for 
central Au + Au (data points) and for p+p plus a model of the Au + Au background (histogram) 
j287j . Note the rotation of the horizontal axis relative to Fig. |2Hl left panel. 



from jet fragmentation in + e~ annihilation, Fig. IHUl shows that the inclusive hadron popula- 
tion in the range of pt*""'^ ~ 4 GeV/c is significantly different in central Au + Au collisions than 
in vacuum jet fragmentation. It therefore cannot be excluded that non-perturbative mechanisms 
contribute to the dihadron phenomenology in Figs. EHlandinZl left panel. Fig. EZl right panel, 
shows the same dihadron correlation analysis described above but for pj'*'''^>6 GeV/c |287j . The 
points are the measured distribution for central Au + Au, in this case not background subtracted, 
while the histogram is the sum of the measured p + p correlation plus a model of the central 
Au + Au background including elliptic flow. The near-angle peaks are again seen to be similar, 
while the back-to-back correlation strength for central Au + Au is again seen to be strongly 
suppressed relative to the p + p excess above background. The persistence of the back-to-back 
suppression to higher px argues that it is indeed due to jet quenching. 

In the direction opposite the triggered hadron, the azimuthal angle distribution measures 
the single inclusive hadron distribution of the away-side jet. The width of the peak in p + p 
collisions is characteristic of the jet profile and is determined by both the distribution of the 
intrinsic pt perpendicular to the jet axis and the relative transverse momentum between the 
two back-to-back jets. In d + Au collisions, initial multiple parton scattering may broaden the 
dijet relative pt and thus the back-to-back dihadron correlation, but the integrated correlation 
strength should remain approximately the same a.s in p + p collisions since the integrated area 
under the away-side peak is essentially determined by the fragmentation function of the jet. 
In central Au + Au collisions, however, partonic energy loss will suppress the leading hadron 
distribution in a jet, thus leading to the reduction of the away-side jet peak as indeed shown by 
the experimental data. 

For the small-angle pairs, the shape of the correlation in p+p collisions is entirely determined 
by the intrinsic pt perpendicular to the jet axis. The integrated value of the correlation strength 
is related to the ratio of dihadron and single hadron fragmentation functions [2EE\. This corre- 
lation remains approximately the same in p + p, d + Au and Au + Au collisions, indicating jet 
fragmentation in all cases. According to the picture of parton energy loss, a parton with reduced 
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energy fragments outside the dense medium, giving rise to a leading dihadron correlation similar 
to that in the absence of energy loss. The lost energy will be carried by the radiated gluons 
which in turn will only contribute to soft hadrons along the direction of the triggered hadron. 

To study the back-to-back dihadron correlation, one can again apply a LO pQCD parton 
model. The spectrum of back-to-back dihadrons from the independent fragmentation of back- 
to-back jets can be calculated as 

E1E2 ^ = IT H / d'^b(frdXadxb(fkaTd'^hTtA{r)tA{\h - r\)gA{kaT, r)gA{hT, |b - r|) 



da 



xfa/A{Xa,Q'^,r)fb/A{Xb,Q'^, |b - r\)Dh/c{Zc, Q"^ , ^Ec)Dh/d{Zd, Q"^ , '^Ed 

i:4 



2Tcz'^zj di 



cd)6 (pa+Pb-Pc-Pd)- (79) 



Let hadron hi be the trigger hadron with pxi = p%^^- We define the hadron-triggered fragmen- 
tation function (FF) as the back-to-back correlation with respect to the triggered hadron: 



D^^^^izT^Jr'^) = ^""1 (80) 



similar to the direct-photon triggered FF in 7-jet events |238[ \'2'A{i\ . Here, zt = Pt/Pt^^ and 
integration over \yi,2\ < is implied. In a simple parton model the dijets will be precisely back- 
to-back, but the initial parton pt distribution will give rise to a Gaussian angular distribution. 
In addition, pt smearing within a jet must be taken into account using a Gaussian distribution 
with a width of ~ 0.6 GeV/c. 

Fig. EHl left panel, shows the calculated back-to-back correlations for charged hadrons in 
Au + Au collisions compared to the STAR data |22Hj . The same energy loss that is used to 
calculate single hadron suppression also describes well the observed back-to-back suppression 
and its centrality dependence. In central Au + Au collisions, multiple parton scatterings that 
induce partonic energy loss can also generate smearing in the transverse momentum of the final 
parton before fragmentation. This can further suppress the back-side correlation at its peak 
|289j . However, after integration over the azimuthal angle, the smearing due to final state 
Pt broadening does not influence significantly the integrated suppression factor of the hadron- 
triggered FF. 

The hadron-triggered fragmentation function is obtained by integrating over 0, D^^^'^ (zt, Pt^^) 
!^/^d(t)D'^^^^{zT,(p^PT'^)- The dihadron suppression factor defined by STAR |223I (Eq. [^) is 



T I trigN _ ^AA V-TiPt ) /g.N 

Iaa{zt,Pt) = T:r777 wiT' (81) 



Jrig _ 



which is the modification factor of the hadron-triggered FF. Fig. EHl right panel, shows the 
suppression factors of the hadron-triggered FF's for different values of p^'^ in central Au + Au 
collisions, compared to a STAR data point obtained by integrating the observed correlation over 
7r/2 < \A(f)\ < IT. The dashed lines illustrate the small suppression of back-to-back correlations 
due to the initial nuclear broadening in d + A collisions. The strong QCD scale dependence of 
the fragmentation functions on p^'^ is to a large extent canceled in the suppression factor. The 
approximately universal shape reflects the weak p^ dependence of the hadron spectra suppression 
factor in Figs. EH due to the unique energy dependence of parton energy loss. 
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Figure 38: Left: Back-to-back dihadron correlations for Au + Au (lower curves) and p + p (upper 
curves), with 4<p^'^<6 GeV/c, 2<pT<p%^^ GeV/c, and \y\ < 0.7, compared to background- 
subtracted STAR data |223j . Right: The suppression factor for hadron-triggered fragmentation 
functions in central (0-5%) Au + Au (d-l-Au) collisions, compared to STAR data |223j . 



Preliminary studies of back-to-back dihadron distributions with low associated px have re- 
vealed a very broad azimuthal distribution at low zt for central Au + Au collisions, with no 
evidence of jet-like azimuthal correlations |290| I291j . The measured distributions are consis- 
tent with statistical momentum balance by a large ensemble of recoiling hadrons even for 
p^^^>6.5 GeV/c |290j . These preliminary results suggest that the soft hadrons in the recoiling 
jet may be significantly broadened, resulting from thermalization of the emitted gluons and pt 
broadening of the jet in the medium. More extensive studies of the modification of the FF at 
low Zt with higher p^'^ are needed to elucidate the picture. 

The dihadron correlation on the same side also places important constraints on coalescence 
models. The same-side correlations characteristic of jet structure exclude the recombination of 
uncorrelated partons from the thermal medium as the dominant mechanism for hadron produc- 
tion at intermediate pt- Such a jet structure on the near side rather favors coalescence of a fast 
parton from the hard scattering with a slower parton from the medium. In the case of hard- 
soft coalescence, both the shape and strength of the correlation may be different from vacuum 
fragmentation. In particular, if the thermal medium develops longitudinal flow, coalescence of 
a leading parton with a thermal parton will lead to broadening of the correlation in rapidity. 



6.7 High Pt Azimuthal Anisotropy 

Non-central collisions generate an initially anisotropic reaction zone (Fig. |22l left panel), with 
the long axis perpendicular to the reaction plane. In pQCD, the final state partons following 
a hard scattering have a priori no correlation with the azimuthal orientation of the reaction 
plane. However, the average path length of parton propagation in the medium will vary with 
the azimuthal angle relative to the reaction plane, leading to an azimuthal dependence of the 
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Figure 39: Left: azimuthal anisotropy V2 from 4-particle cumulant analysis of 200 GeV Au + Au 
collisions from STAR |i297j compared to parton model calculations incorporating partonic energy 
loss. Right: pt dependence of V2 for identified particles, with both axes scaled by the number 
of constituent quarks; data from PHENIX and STAR, figure from |298| . 



total partonic energy loss. Thus, azimuthal anisotropy or "elliptic flow" of high hadrons 
relative to the reaction plane is a consequence of partonic energy loss and is insensitive to effects 
of initial state interactions |293t I236j . 

The experimental techniques for studying azimuthal anisotropy at high pj- are the same 
as at low Pt (Sect. 15. 3|) . The strength of the correlation is measured by the second Fourier 
coefficient V2 of the azimuthal distribution [Eq. ()41|) ]. Special care must be taken to account for 
so-called "non-flow" effects contributing to V2 which arise from multiparticle correlations that are 
unrelated to correlations with the reaction plane |158j . Non-flow effects may result, for instance, 
from resonance decays, momentum conservation, and particularly at high px, intra- and inter-jet 
hadron correlations. The standard method of correlating each high pt hadron with the reaction 
plane is equivalent to a dihadron correlation analysis [158:. Higher order cumulants have been 
shown to be markedly less sensitive to non-flow effects than the two-particle correlation methods 
|2941 1295j . It has recently been argued, however, that higher order cumulants are susceptible 
to fluctuation effects and may thereby underestimate the true correlation strength |29tij . The 
two-particle and higher order correlation measurements thus may bracket the true flow V2 and 
are usually reported together. 

Fig. I2H1 right panel, shows the azimuthal anisotropy f 2 in non-central Au + Au collisions 
at 200 GeV for the reaction plane (circles) and 2- (triangles) and 4-particle (stars) cumulant 
methods |224j . The measured V2 is large, in the sense that its magnitude is similar to the 
eccentricity due to the initial spatial anisotropy of the collision [Eq. (jH^ ] |299| 1272] . Since the 
spatial eccentricity of the flreball is diluted by expansion |113t 1139] , such a large f 2 at high pt 
may result from strong partonic energy loss at the earliest, hot and dense phase of the evolution. 

In the parton model approach, the azimuthal dependence of the partonic energy loss for 
non-central heavy- ion collisions is given by Eq. ()76|) . The anisotropic energy loss in the effec- 
tive modifled fragmentation functions of the pQCD parton model can be used to obtain the 
azimuthally anisotropic hadron spectra at high px- Fig. [HHl left panel, shows V2 for charged 
hadrons generated from partonic energy loss (dot-dashed) compared to the measured f 2 from a 
4-particle cumulant analysis |224t I297j . The energy loss extracted from high pt inclusive sup- 
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pression can account for the observed azimuthal anisotropy at large pt>Q GeV/c (dot-dashed 
hne) . 

The azimuthal anisotropy at high pt can also depend on the transverse expansion, which 
depletes the gluon density more rapidly than one-dimensional expansion and dilutes the ini- 
tial geometric anisotropy. The depletion in gluon density is however compensated by longer 
propagation time in the medium, giving rise to roughly the same total energy loss. In realistic 
calculations the azimuthal anisotropy of partonic energy loss is found to decrease slightly rela- 
tive to the case of no transverse expansion |26()j and quantitative analyses should take this effect 
into account. Note also that in the current parton model calculations a hard-sphere nuclear 
distribution is used. A more accurate calculation would employ the more realistic Wood-Saxon 
distribution which would reduce the high pt V2 |3UUj . However, partonic energy loss occurs 
mainly at early time due to rapid expansion, making the final result relatively insensitive to the 
shape of edge of the reaction zone. 

The observed f 2 at intermediate pj- is larger than the simple parton model calculation. This 
discrepancy may be attributable to effects of parton coalescence. If the difference is due to 
the flow of kaons and baryons generated by coalescence, they must have V2 ~ 0.23 for 20-50% 
collisions and V2 ~ 0.11 for 0-10% collisions. The total V2 incorporating this effect is shown by 
the solid lines in Fig. EHl left panel. 

A remarkable phenomenological scaling for particle-identified f 2 is seen in Fig. EHl right 
panel, where both the horizontal (pt) and vertical (^2) axes have been scaled by the number of 
constituent quarks in the hadron {ug = 2 for mesons, Ug = 3 for baryons). The scaled distribu- 
tions for all hadrons except pions collapse within experimental uncertainties to a universal curve, 
indicating that constituent quark dynamics may be driving the V2 of hadrons at intermediate 
Pt- Such scaling emerges from a recombination model based an exponential partonic spectrum 
|30H I26[J] . with the requirement that all partons contribute equally to the hadron's momentum 
(i.e. the Wigner function of the hadron is a (5-function in momentum space). Relaxation of 
the 5-function condition is found however to generate negligible change in V2 at intermediate 
P'r |30H 1269] ■ so that the observed scaling provides robust support for recombination as the 
mechanism underlying hadronization at intermediate pt. Pions are excused from following the 
scaling because they arise dominantly from p-decay |302j . 

If hadrons at intermediate pt are produced through recombination of a fast parton from a 
jet and with soft partons from the thermal medium, their azimuthal anisotropy relative to the 
reaction plane should depend both on the azimuthal dependence of the partonic energy loss and 
on the elliptic flow of the soft partons. In such models the functional dependence of f 2 on p^ for 
soft partons must be specified, with parameters fixed by fitting to the data. The predictive power 
of the model lies in the flavor dependence. Since baryons receive a larger relative contribution 
from coalescence, with the effect extending to larger pt than for mesons, splitting of V2{pt) for 
baryons and mesons should be observed. Fig. EUl shows measured V2{pt) for K and A compared 
to two coalescence model calculations |3011 1276] . The models reproduce f 2 of all species well up 
to pr~4-5 GeV/c. As for inclusive spectra, coalescence effects on V2 diminish for px > 5 GeV 
where V2{pt) for all hadron species should become the same, driven by partonic energy loss. 

An additional test of partonic energy loss results from the differential study of high pt 
dihadron correlations relative to the reaction plane orientation. Fig. ^2 shows a preliminary 
study of the high p^ dihadron correlation for non-central (20-60%) Au + Au collisions, with 
the trigger hadron situated alternatively in the azimuthal quadrants centered on the reaction 
plane ( "in-plane" ) or those orthogonal to it ( "out-of-plane" ) |303| I304j . The same-side dihadron 
correlation in both cases is similar to that in p + p collisions. In contrast, the suppression of the 
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Figure 40: V2 of K and A |146j compared to coalescence model calculations. Left: calculations 
from |276j . Model predictions for and Q also shown. Dotted curve is strange quark anisotropy. 
Right: calculations from |301j . 

back-to-back correlation depends strongly on the relative angle between the trigger hadron and 
the reaction plane. This systematic dependence is also consistent with the picture of partonic 
energy loss: the path length in medium for a dijet oriented out of the reaction plane is longer 
than in the reaction plane (Fig. |221 left panel), leading to correspondingly larger energy loss. 

6.8 Partonic Energy Loss vs. Hadronic Absorption 

From the analyses within a parton model of single inclusive and dihadron spectra and the 
azimuthal anisotropy V2{pt), the average energy loss for a 10 GeV quark propagating through 
the expanding medium created at initial time Tq = 0.2 fm/c in 200 GeV central Au-|-Au collisions 
is {dE/dL)id ~ 0.85 ± 0.24 GeV/fm. This is equivalent to dEo/dL ^ 13.8 ± 3.9 GeV/fm in a 
static and uniform medium over a distance Ra = 6.5 fm at initial time Tq = 0.2 fm. The value is 
about 30 times higher than the quark energy loss in cold nuclei, as extracted from HERMES DIS 
data. Since the partonic energy loss in the thin plasma limit is proportional to the gluon number 
density, this indicates that the initial gluon density at tq = 0.2 fm/c reached in central Au + Au 
collisions at 200 GeV is about 30 times higher than the gluon density in a cold Au nucleus. This 
number is consistent with the estimate from the measured rapidity density of charged hadrons 
|l()lj using the Bjorken scenario jll2j . assuming duality between the number of initial gluons 
and final charged hadrons. Since total parton energy loss is only sensitive to the gluon density 
of the medium [Eq. (jU^ ]. this is the only property that can be extracted from the suppression 
of single and dihadron spectra. To extract the initial energy density, additional measurements 
such as jet broadening are required. Alternatively, the energy density can be estimated from 
global measurements. Given the measured total transverse energy dET/drj ^ 540 GeV or about 
0.8 GeV per charged hadron in central Au + Au collisions at ^/s = 130 GeV |134j . the initial 
energy density is 50-100 times that in cold nuclear matter (Sec. 14. 5|) . 

In the above analyses of RHIC data, the mechanism dominantly responsible for jet quenching 
is partonic energy loss prior to hadronization of the jet. While this picture is in good accord 
with the observed high pj- phenomena, it is essential to ask whether hadronic interactions, 
specifically the interaction of hadronic jet fragments with the medium, can at least in part 
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Figure 41: Background-subtracted high px dihadron correlation for different orientations of the 
trigger hadron relative to the reaction plane ^30,^ I3()4j . 

generate the observed high pt phenomena and contribute substantially to the jet quenching 
|3U5| I3U61 l3U7j . Some simple considerations already argue against this scenario. The formation 
time of hadrons with energy Eh and mass ruh is = {Eh/mh)Tf, where the rest frame formation 
time Tf ~ 0.5 — 0.8 fm/c. Thus, a 10 GeV/c pion has formation time ~ 50 fm/c and is unlikely 
to interact as a fully formed pion in the medium. Since the formation time depends on the 
boost, the suppression due to hadronic absorption with constant or slowly varying cross section 
should turn off with rising p^, at variance with observations (Fig. |^ right panel). A detailed 
hadronic transport calculation j307j leads to a similar conclusion: the absorption of formed 
hadrons in the medium cannot account by a large factor for the observed suppression, and the 
suppression is attributed to medium interactions of "pre-hadrons" which have short formation 
time and constant cross section, in other words properties similar to those of colored partons 
|307j . Further consideration of the available high px data j308j also supports the conclusion that 
jet quenching in heavy ion collisions at RHIC is the consequence of partonic energy loss. In 
particular, large V2 at high px and the systematics of the small-angle dihadron correlations are 
difficult to reconcile with the hadronic absorption scenario. 

Azimuthal anisotropy of the spectra in non-central collisions arises from the initial spatial 
eccentricity of the dense medium. The eccentricity decreases rapidly with time due to transverse 
expansion |113l I131J] . so that momentum anisotropy must be generated early, ~few fm/c. The 
finite anisotropy observed at high px (Fig. |2Hl right panel) is therefore unlikely to have significant 
contribution from the absorption of formed hadrons. 

As we have shown in Fig. |2Hl left panel, the back-to-back dihadron correlation is suppressed 
in central Au + Au collisions relative to p + p and d + Au collisions, while the same-side cor- 
relations are similar in all three systems. In the framework of jet fragmentation, the same-side 
correlation measures the conditional distribution of the second leading hadron within a jet in 
coincidence with the trigger hadron, essentially the ratio of dihadron to single hadron fragmen- 
tation functions. Insofar as both hadrons originate from the same fragmenting parton, this ratio 
will largely be independent of the parton energy. Gluons radiated via energy loss will generate 
soft hadrons and will enhance the conditional distribution at small z. A recent study of the 
dihadron fragmentation functions shows that the conditional (or triggered) dihadron distribu- 
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Figure 42: Average transverse energy [Et)^^^ of the initial partons that produce a final hadron 
with function of (A^part) for different values of p^^^ (lowest curves for smallest px)- 



tion at large z within a jet is quite stable against radiative evolution (or energy loss), while soft 
secondary hadrons are enhanced f288j . 

Hadronic absorption, on the other hand, will generally suppress the leading and secondary 
hadron independently, thereby suppressing the conditional distribution to the same degree as 
the single particle inclusive yield. This may not be the case if the hadronic absorption is very 
strong and only jets originating in a thin surface shell generate trigger hadrons which are not 
suppressed, but we consider this scenario to be unrealistic. Calculations in a parton model with 
energy loss show that even the trigger-biased jet population generated near the surface loses on 
average about 2 GeV of energy j208j, as seen in Fig. |121 This energy will be carried by soft 
hadrons correlated with the trigger, and the total jet energy for fixed trigger hadron p^'^ should 
be larger in Au + Au collisions than m p + p collisions by about 2 GeV. We will discuss the 
prospects for observing the radiated energy carried by soft hadrons using dihadron measurements 
in section inini 

If hadronic absorption suppresses high px hadrons and jets, it should also do so in heavy-ion 
collisions at SPS energy. Hadronic spectra at this energy vary strongly with px and are very 
sensitive to initial transverse momentum broadening and parton energy loss [HH] . The measured 
7r° spectrum in central Ph + Ph collisions appears to exhibit only the expected Cronin enhance- 
ment relative to a p + p reference, with no suppression observed |30^l31Ull311j . However, the vr" 
yield for central collisions is seen to be suppressed relative to that for peripheral collisions |31Uj . 
It has been been pointed out recently |312j that the p + p reference used in this study contains 
significant uncertainties, and a reassessment also reveals a possible high px vr° suppression for 
central Ph + Ph relative io p + p collisions. Recent analysis of dihadron correlations shows that 
both same-side and back-to-back jet-like correlations are not suppressed in central collisions at 
the SPS, though the back-side distribution is broadened |313j . The question of high px hadron 
suppression at the SPS and its connection to jet quenching therefore remains open. More gen- 
erally, study of the dependence of the suppression will provide an essential cross-check of 
our understanding of these phenomena. High px hadron production results from the recently 
completed 62 GeV Au + Au run at RHIC are now becoming available |314j . 
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We conclude that the data provide no clear support for hadronic absorption as the dominant 
mechanism underlying the observed high pt suppression phenomena. We consider large partonic 
energy loss in the dense medium formed in central Au + Au collisions at RHIC to be well 
established. We now discuss future measurements that will exploit this new phenomenon as a 
precision probe of the medium. 

6.9 Hard probes: an outlook 

The discovery of jet quenching opens a new era in the study of dense QCD matter. Higher 
precision data extending to larger p^ will map in detail the modification of the fragmentation 
functions and the energy dependence of the energy loss. The ultimate measurement of modified 
FF's in heavy ion collisions will be the direct photon triggered FF. Measurements of charmed 
mesons will also provide key tests of the parton energy loss scenario of jet quenching; recent 
theoretical studies reveal features of heavy quark energy loss that are measurably different from 
those of light quarks and gluons j315| I316| I317| I318j . Dihadron fragmentation and correlations 
within the same jet will also be important, for instance the broadening of the dihadron angular 
correlation which probes the average momentum transfer of the interaction with the medium 
that may be related directly to the energy density |288j . 

Direct Photon Tagged Jet Quenching 

The definition of direct photon-triggered fragmentation functions is similar to the hadron- 
triggered fragmentation function in Eq. ()80p . replacing the triggered hadron with a direct photon 
|238| I239j . This corresponds to the measurement of the hadron pr distribution in the opposite 
azimuthal direction to the triggered photon. Since a direct photon in the central rapidity region 
{y = 0) is always accompanied by a jet in the opposite azimuth with roughly equal transverse 
energy, the pt distribution of particles is directly related to the fragmentation function of a jet 
with known initial energy, E"^* pa E^. By comparing the extracted jet fragmentation function 
in A + A to that in p + p collisions, the modification of the fragmentation function and thereby 
the partonic energy loss can be measured directly. 

Fig. 1^ shows examples of the the suppression factors of the direct photon-tagged jet fragmen- 
tation function for average partonic energy loss dE/dx = 1 GeV/fm. They are similar in shape 
to the hadron-triggered fragmentation function. The presently measured single particle inclusive 
and dihadron suppression generate an average energy loss of about dE/dx = 0.85 GeV/fm for 
a 10 GeV quark in central Au + Au collisions at RHIC. According to the parameterized energy 
dependence in Eq. (jUHj) . the average energy loss for a 15 GeV quark should be about 1.26 GeV, 
slightly larger than that used to generate the curves in the figure. The initial mean free path 
that fits the inclusive and dihadron data, Aq = 0.3 fm, is consistent with 1-2 fm used in the 
calculation after correction for the one-dimensional expansion. 

The measurement of direct photon-triggered fragmentation is extremely challenging because 
of low cross sections and the difficulty of isolating direct photons in the heavy-ion collision 
environment. However, statistical measurements of the inclusive direct photon yield are now 
becoming available. Fig. H^l right panel, shows the preliminary measured ratio of the yield of 
total (direct + decay) photons to expected decay photons, from PHENIX j319j . The measured 
ratio in p+p collisions is consistent with an NLO calculation of the direct photon yield |319j . The 
figure compares the ratio in central Au + Au collisions to NLO calculations of the direct photon 
yield normalized alternatively by a background derived from the binary-scaled 7r° yield from p+p 
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Figure 43: Left: Calculated modification factor of the photon-tagged jet fragmentation function 
in central Au + Au collisions at ^/s = 200 GeV for fixed dEg/dx = 1 GeV/fm [239J. Right: Ratio 
of total to decay photon yield in central Au + Au collisions at 200 GeV, from PHENIX |319j . The 
lines correspond to direct photon yields from NLO pQCD and decay photons from the measured 
7r° spectra. Upper curves incorporate Au + Au jet quenching in the background estimate, lower 
curves are for binary-scaled background. 



collisions or by the measured, highly suppressed vr yield from central Au + Au collisions. The 
latter normalization is clearly favored, demonstrating unambiguously that the photon excess 
is not suppressed. This is consistent with the partonic energy loss picture since the photon 
does not carry color charge. Further isolation of direct photons from the QCD fragmentation 
background is difficult in heavy ion collisions because standard techniques, in particular isolation 
cuts, cannot be directly applied. Strong jet suppression may help in this regard, however. This 
result in any case shows great promise and measurement of the modified fragmentation function 
of a jet recoiling from a hard photon is on the horizon. 

Heavy Quark Energy Loss 

The current experimental studies of jet quenching have focused on light quark and gluon jets. 
The energy loss of gluons and light quarks differ by a factor of Ca/Cf = 9/4, but it is difficult to 
differentiate quark and gluon jets in practice in p + p and Au + Au collisions. Identified hadron 
ratios such as K~ / or p/p may be used to tag contributions from different flavor jets |320tl32lj 
but this does not provide sharp discrimination, due to flavor mixing in the fragmentation of 
light quark jets. The study of heavy quark jets offers cleaner flavor tagging. Recent studies have 
revealed observable features of heavy quark energy loss in medium that may provide signiflcant 
new insight into the partonic energy loss mechanism |315| 13161 1317| I318j . 

The formation time of gluon radiation off a heavy quark, measured with respect to the 
propagation of the heavy quark inside the medium, is reduced relative to a light quark because 
of the large quark mass. LPM interference should therefore be reduced signiflcantly for heavy 
quarks at intermediate energy. In addition, the heavy quark mass suppresses the gluon radiation 
amplitude at angles smaller than the ratio of the quark mass to its energy j315j . Both mass 
effects will lead to different energy loss for a heavy relative to a light quark in a dense medium. 

The mass dependence of the gluon radiation amplitude comes from the heavy quark prop- 
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agators. It suppresses the induced gluon spectrum for small angle radiation for a heavy quark 
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Here 6 = It/q~z and 9q = M/q~ represent the radiation angle and the "dead-cone" within 
which the gluon radiation is suppressed. This suppression leads to a reduced radiative energy 
loss of a heavy quark. To illustrate the mass suppression of radiative energy loss imposed by 
the "dead-cone", the ratio {Az^){xb, Q"^)/ {^z'^){xb, Q^) of charm quark and light quark energy 
loss in DIS off a nucleus is plotted in Fig. EH left panel, as functions of xb- 

Within the same framework of DIS off nuclei, the gluon formation time for radiation from a 
heavy quark 

^~P+xl 4 + (1-z)2M2' ^ ' 

is shorter than that for gluon radiation from a light quark, with significant consequences for 
heavy quark energy loss. Because of the quark mass dependence of the formation time relative 
to the nuclear size in the given frame, mi^RA/rfp^ ~ x bM"^ / x aQ"^ ■, there are two distinct 
limiting behaviors of the energy loss for different values of xbIQ'^ relative to xa/M'^. When 
xb/Q"^ ^ xa/M'^ for small quark energy (large x^) or small Q^, the formation time of gluon 
radiation off a heavy quark is smaller than the nuclear size. In this case, there is no destructive 
LPM interference. The heavy quark energy loss 

is linear in nuclear size Ra- In the opposite limit, xb/Q"^ ^ xa/M"^, for large quark energy 
(small Xb) or large Q^, the quark mass becomes negligible. The gluon formation time could still 
be much larger than the nuclear size. The LPM interference will limit the available phase space 
for gluon radiation and the heavy quark energy loss 

now has a quadratic dependence on the nuclear size, similar to the light quark energy loss. 
Fig. right panel, shows numerical results for the Ra dependence of charm quark energy 
loss rescaled by C{Q'^)CA0i1{Q'^)/Nc, for different values of xb and Q^. The Ra dependence is 
quadratic for large values of or small Xb-, but becomes almost linear for small or large xb- 
The charm quark mass is set at M = 1.5 GeV in the numerical calculation. 

Apparently, energy loss induced by gluon radiation is significantly suppressed for heavy 
relative to light quarks when the momentum scale Q or the quark initial energy q~ is not large 
compared to the quark mass. Only in the limit M <^ Q, q~ , is the mass effect negligible and the 
heavy quark energy loss approaches that of a light quark. In the kinematic regime accessible at 
RHIC, the heavy quark energy loss will be much reduced relative to the light quarks and gluons. 
In addition, the mass effect is likely to reduce the effect of thermal absorption of gluons, whose 
average energy is much smaller than the heavy quark mass. 

Heavy quark physics at top RHIC energy and design luminosity is studied primarily through 
charm production, due to the relatively low rate of beauty quark production. Several techniques 
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Figure 44: Left: The xb dependence of the ratio between charm quark and hght quark energy 
loss in DIS off a heavy nucleus. Right: Dependence of charm quark energy loss on nuclear size 
Ra, for different values of and xb- 



are available to study charm production directly and indirectly. Fig. |33 left panel, shows the 
spectrum of electrons from the semi-leptonic decay of D mesons in 130 GeV Au + Au collisions 
from PHENIX [322; , which provides an indirect measurement of the open charm meson spectrum. 
Preliminary non-photonic electron spectra from 200 GeV p + p, d + Au and Au + Au collisions 
have also been reported [323] . Within the current uncertainties the charm yield in all three 
systems appears to scale as the number of binary collisions, though the data are also consistent 
with small suppression at large px in Au + Au. 

Single electron spectra resulting from semi-leptonic decays are however rather insensitive to 
changes in the D-meson spectrum |324j and more precise measurements of charm quark energy 
loss require direct measurement of fully reconstructed D decays into hadrons. STAR recently 
reported direct reconstruction of D-mesons in d + Au collisions. Fig. 021 right panel, shows a 
preliminary D meson spectrum constructed from various hadronic decay channels |325j . The 
spectrum is consistent with the single electron spectrum measured by the same experiment. Such 
measurements are considerably more demanding in Au + Au collisions but will provide direct 
study of the open charm production and the effect of heavy quark energy loss. 

Parton recombination may also play a role in charm production. Such effects have already 
been shown to be significant in the forward direction for h+p collisions |267j . Since fragmentation 
functions are harder for heavy than light quarks, heavy quarks carry a large fraction of the heavy 
meson momentum in the wave-function, making it easier for a fast heavy quark (anti-quark) to 
pick up a slow light anti-quark (quark) to form a heavy meson. The recombination effect is 
therefore expected to be more significant for D mesons than light quark hadrons. As in the 
intermediate pt region of light hadron spectra, discrimination of the effects of partonic energy 
loss and parton recombination may be accomplished most clearly via azimuthal anisotropy of 
the D meson spectra. Some recombination models have predicted small but finite azimuthal 
anisotropy of the D meson spectra due to recombination with thermal light quarks |326j . This 
measurement is however very demanding and most likely requires upgrades to the present RHIC 
detectors. 
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Figure 45: Left: Non-photonic single electron spectra in Au + Au collisions at 130 GeV, from 
PHENIX |322j . Right: Spectra of fully reconstructed D meson spectra in ci + Au collisions at 
200 GeV from STAR [2H^ . 



Dihadron Fragmentation Functions and Angular Correlations 

The suppression of single inclusive hadron spectra or the modification of the single hadron 
fragmentation function provide direct measurement of the gluon density of the dense medium. 
Broadening of the jet cone in principle could probe the average transverse momentum transfer to 
the jet parton, which can be used to directly estimate the energy density rather than the gluon 
density of the medium |327| 1328] . Such a study may be possible using dihadron correlations. 
Dihadrons can result from the hadronization of a single leading parton or from two independent 
partons, one leading and the other a radiated gluon. The dihadron fragmentation or correlation 
function can therefore probe multiple scattering and induced gluon radiation in the medium. The 
first mechanism dominates when both hadrons carry a large fraction of the initial parton energy. 
This is likely the case for the same-side correlations measured by STAR fFig.l^Hj) and may explain 
why the correlation does not change despite the energy loss expected to be suffered by the leading 
parton. If the associated (secondary) hadron carries smaller fractional momentum it could come 
from radiated gluons. In this case the dihadron fragmentation function should be enhanced, as 
seen by a recent study of its evolution ^288^. Due to the transverse momentum transfer from 
multiple scattering, the angular correlation between the hadrons should be also broadened. This 
is best seen by looking at the soft hadrons in a jet, though precision measurements of such 
observables are extremely challenging. 

A first attempt to measure the radiated energy from the associated hadrons (Fig. is 
seen in Fig. EHl left panel, which shows the integrated multiplicity and summed "scalar pt" 
of hadrons within 0.15 < < 4 GeV/c correlated with a trigger with pt*"^>4 GeV/c |291j . 
The integration is performed over phase space broader than the normal jet size in vaccum, with 
"near-side" corresponding to |A</)|<1.1 and |A?7|<1.4 and "away-side" corresponding to |A0|>1.1 
and |?7|<1.1. Background yield is subtracted assuming that the correlated yield is negligible in 
the region 0.9<| A(/)|<1.3. The associated multiplicity and summed pt above background show 
increases from p + p to Au + Au collisions qualitatively similar to the theoretical expectations 
from partonic energy loss, though the moderate trigger p^ makes direct comparison with such 
calculations difficult. Future analyses will increase the trigger pt and reduce the systematic 
uncertainties associated with the background definition. Effects of parton recombination in the 
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I It (GeV/c) 

Figure 46: Left: charged hadron multiplicity (upper) and total scalar px (lower) of the near and 
away-side dihadron distributions relative to a trigger hadron, integrated over broad (Ar^, A0) 
intervals for 4 < p^*^ < 6 GeV/c and 0.15 < pt < 4 GeV/c in p + p (open symbols) and Au + Au 
collisions, from STAR |291j . See text for details. Right: Charged hadron Jt distributions [Eq. 
(jSni)] from fully reconstructed jets in p+p coUisions at ^/s=200 GeV from STAR |H29j . compared 
to simulations based on PYTHIA . The two distributions have different lower bounds 

on the hadron pt relative to the beam axis. 
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jet structure must also be addressed in order to extract net effect of partonic energy loss. 
Pull jet reconstruction 

As we have discussed extensively, the study of correlations among jet fragments on a statistical 
basis provides key observables for studying hard probes in collisions of heavy nuclei. Event-wise 
full jet reconstruction with good energy resolution is exceedingly difficult in all but perhaps the 
most peripheral collisions of heavy nuclei, due to the large multiplicity and the complexity of the 
underlying event. Standard jet reconstruction techniques can however be applied to collisions of 
simpler systems such as p + p and d + Au. Fig. 0^1 right panel, shows the preliminary results 
on the Jt distribution in fully reconstructed jets in 200 GeV p + p collisions j329j . where Jt is 
the component of hadron momentum perpendicular to the jet thrust axis: 



Ph is the hadron momentum and Pjet is the jet momentum. The jets are at mid-rapidity and have 
measured transverse energy (Ej-) ~ 11 GeV. The figure shows comparison to PYTHIA-based 
simulations |330j which describe the data well, including variations in the Jt distribution due to 
variations of the lower limit of the hadron pt relative to the beam axis (the kinematic "seagull 
effect"). Aside from the importance of jets for studying QCD processes at RHIC energies, jet 
measurements in simpler systems at RHIC (potentially including the collisions of light nuclei) 
provide essential data-based calibrations of the more limited multi-hadron correlation observables 
that are accessible in the collisions of heavy nuclei |H29j . 

7 Summary 

The Relativistic Heavy Ion Collider has initiated a new era in the study of QCD matter under 
extreme conditions. The RHIC accelerator has unprecedented flexibility for a collider, having 
generated significant integrated luminosity for Au + Au collisions at several energies as well 
as for polarized protons and the very important d + Au control experiment. The four RHIC 
experiments have produced a large body of high quality data. There is considerable overlap in 
their physics coverage and it is notable that the results are in agreement on all major physics 
points. 

The data collected and analysed in the first three years of RHIC operations indicate that a 
dense, equilibrated system is generated briefly in the most violent, head-on collisions of heavy 
nuclei at top RHIC energies. The most economical explanation of the observed phenomena is that 
a state of matter dominated by colored (partonic) degrees of freedom has been produced, though 
direct observation of the deconfinement transition awaits further experimental and theoretical 
developments. For long wavelength excitations (low or soft probes) the matter evidently 
responds as a near-ideal, strongly coupled fluid, while for short wavelength (high Q"^) probes 
it is highly dissipative. The medium is markedly different from the ideal non-interacting gas 
expected from QCD at asymptotically high temperature. 

In this review we have presented only a limited subset of the available results, concentrating 
on the most mature measurements that directly address new phenomena. The evidence that a 
novel state of matter has been created is based on several lines of argument (see also I332j ): 




(86) 
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• Collective behaviour: the broad success of hydro dynamic calculations in describing the 
inclusive spectra and azimuthal anisotropies of soft hadrons indicates that local thermal 
equilibrium is established at time t<1 fm/c after the collision and that the matter expands 
as a near-ideal fluid. 

• Partonic energy loss: the strong suppression of inclusive yields and correlations of high 
Pt hadrons shows that high energy partons dissipate significant energy in the medium, 
providing direct evidence of strong interaction among partons that is essential to establish 
thermalization. The absence of high suppression in the d + Au control experiment 
confirms that the STipprcssion is due to final state interactions in dense matter. Hadronic 
absorption calculations cannot account qualitatively for the systematic behavior of the 
phenomena, and the suppression agrees with a picture of partonic interactions with a 
medium having high color charge density. 

• Energy density: the independent analyses of transverse energy production, hydrodynamic 
flow, and partonic energy loss result in a consistent estimate of the energy density early in 
the collision of 5-10 GeV/fm^, well beyond the critical value for transition to a deconfined 
phase expected from Lattice QCD. 

Nevertheless, there remain important open questions. The mechanisms underlying the ap- 
parent very rapid thermalization arc at present unclear. The role of finite viscosity, especially 
at the later hadronic stage of the expansion, remains to be fully understood. Most importantly, 
sensitivity to the Equation of State has not yet been systematically explored. New data are crit- 
ical for addressing these issues, in particular more precise measurements of the low pr spectra 
and asymmetries of multistrange baryons and charmed mesons. New theoretical developments 
and systematic studies that are constrained by the broad range of data now available are likewise 
needed. 

In the area of hard probes, striking effects have been observed but precision data in a larger 
Pt range are needed to explore the detailed properties of jet quenching and their connection to 
other properties of the dense matter. The region 2<pt<6 GeV/c has significant contributions 
from non-perturbative processes, perhaps revealing novel hadronization mechanisms. However, 
all studies to date of azimuthal anisotropies and correlations of "jets" have by necessity been 
constrained to this region, with only the inclusive spectra extending well beyond the range where 
non-perturbative processes are seen to play a role. High statistics data sets for much higher px 
hadrons are needed to fully exploit azimuthal asymmetries and correlations as measurements of 
partonic energy loss. Heavy quark suppression is theoretically well controlled, and measurement 
of it will provide a critical check on our understanding of partonic energy loss. The differential 
measurement of energy loss through measurement of the emerging away-side jet and the recovery 
of the energy radiated in soft hadrons is still in its initial phase of study. A complete mapping 
of the modified fragmentation with larger initial jet energy and with a direct photon trigger 
will cross check the energy dependence of energy loss extracted from single inclusive hadron 
suppression. Experiments at different colliding energies are also essential to reveal the onset of 
critical phenomena, or at the minimum to map the variation of jet quenching with initial energy 
density and the hfetime of the dense system. 

Qualitatively new observables are also on the horizon. With accumulated luminosity in- 
creasing each year, measurements of J/ip suppression in Au + Au collisions will become more 
quantitative and differential. Together with d + Au results and the open charm measurements it 
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may be possible to disentangle final state suppression, initial state nuclear absorption, and pos- 
sibly the contribution to J / ip production from charm quark recombination in the dense medium 
\6'6'6\ These measurements will provide an independent measurement of the initial condi- 

tions of the dense medium. Direct thermal photon emission from the plasma at intermediate pt 
could in principle provide direct measurement of the effective temperature or average parton en- 
ergy of the interacting dense matter. Emission from the hadronic phase and photon production 
from perturbative hard processes form large backgrounds, however, making the extraction of the 
thermal photon yield from the QGP phase extremely difficult [UlinnSl- Finally, experimental 
study of dilepton spectra in the low and intermediate mass region will provide vital information 
on the medium modification of vector mesons due to chiral symmetry restoration, as well as the 
thermal dilepton emission from the plasma phase. 

There has been considerable recent interest in universal properties of QCD at very low 
Bjorken xsj Due to the non-Abelian nature of QCD, the gluon density in hadrons grows 
rapidly with decreasing xbj- However, the gluon self-interaction must eventually lead the gluon 
density to saturate. In the saturation region the system is dense but weakly coupled, controlled 
by the large saturation scale. This is a unique QCD regime, called the Color Glass Condensate 
(CGC). While arguments supporting the CGC are generic, the momentum scale at which it 
occurs must be fixed by measurements. CGC effects are expected to be amplified in heavy nuclei 
and evidence for them has been sought in RHIC data. The question at present is still open, but 
it is possible that the CGC is the underlying initial coherent state of partons in the nucleus from 
which the incoherent, thermalized QGP emerges following the collision. It has been suggested 
|262j that measurements in the forward {d direction) region of c? -|- Au collisions could provide a 
crucial test of saturation phenomena, in particular its spill-over to the large pt region, though 
that region is known to be susceptible to other non-perturbative effects. Measurements in the 
central rapidity region at LHC energies may provide clearer evidence of saturation phenomena. 
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